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INTRODUCTION 



The physics of folding encompasses a wide variety of topics ranging from biology to 
pure mathematics. Folding problems arise for instance in the study of protein confor- 
mations 0] and in the physical modelling of biological or artificial membranes 0. In a 
different context, the physics of paper folding also displays interesting phenomena, in con- 
nection with the mechanical and statistical properties of thin elastic plates or membranes 
0. At a more abstract level, the very description of the folding degrees of freedom of a 
surface gives rise to very interesting statistical models and quite involved combinatorial 
problems. 

In this wide subject, we choose to concentrate on this last type of questions by focusing 
on statistical models of discrete lattice folding and emphasizing the crucial role played by 
geometrical constraints on the statistics of large folded objects. More precisely, we will 
consider discrete models of "solid" or "fluid" membranes in the form of folding problems 
for respectively fixed and random lattices. An example of flxed lattice is a piece of the two- 
dimensional triangular lattice, viewed as a regular collection of rigid equilateral triangles 
whose edges may serve as hinges between neighboring triangles. An example of random 
lattice is a triangulation made of the same triangles arranged into a possibly irregular and 
fluctuating surface tessellation, allowing for curvature defects (more or less than 6 triangles 
around a node) or topological defects (tessellations of a surface with handles). 

Our flrst task will be to characterize the folding degrees of freedom in these models. 
This happens to be an extremely difficult problem in the presence of self-avoidance. For 
this reason, we shall flrst consider the much simpler "phantom" folding problem in which 
the membrane is allowed to interpenetrate itself. A flrst description uses local normal 
vectors as effective spin variables subject to various geometrical constraints. Although 
local, these constraints induce a non-local behavior with propagating creases, resulting in 
a rich phase diagram even for regular lattices. 

A second and more powerful description uses loop conflgurations on the lattice and 
allows to view folding as a particular instance of the more general problem of the "fully 
packed" loop gas, that is a set of non-intersecting loops visiting all nodes of the lattice. 
This approach, combined with fleld-theoretical techniques, allows for a number of predic- 
tions on the statistics of large folded lattices. For regular lattices, these predictions rely 
either on exact solutions via Bethe Ansatz (integrable case) or on effective fleld-theoretical 
descriptions via Coulomb gas. 
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In the case of random lattices, the loop gas interpretation leads in particular to precise 
predictions for various configurational exponents characterizing the asymptotics of the 
number of allowed folding or loop configurations. These predictions make extensive use of 
interpretations of the folding configurations in terms of discrete two-dimensional quantum 
gravity, namely statistical models defined on dynamical tessellations of surfaces. At large 
scales, these models are described by conformal field theories coupled to two-dimensional 
quantum gravity, and the configuration exponents are obtained from a proper identification 
of the central charge of the conformal field theory at hand. Alternatively, the discrete 
models of folding may be represented via diagrammatic expansions of matrix integrals, 
allowing for a number of direct computations. 

Returning to the crucial question of self-avoidance, we will concentrate on the one- 
dimensional self-avoiding folding problem, also known as the "stamp folding problem". 
This problem happens to be equivalent to a famous combinatorial problem, the meander 
problem which may be stated as follows: "enumerate all topologically inequivalent con- 
figurations of a non-intersecting circuit crossing a straight river through a given number 
of bridges". Remarkably, in its meander formulation, the one- dimensional self-avoiding 
folding problem is described by a two-dimensional gas of fully packed loops on a random 
lattice, hence belongs to the class of two-dimensional phantom foldings of random lattices. 
A remarkable outcome of our study is the predicted value for the meander configuration 
exponent 



a = (29 + V145)/12 (1) 
which governs the large n asymptotics of the number of meanders with 2n bridges oc 

This review is based mainly on a series of articles by the authors [fflTTll . Other reviews 
may be found in 
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PART A: FOLDING OF REGULAR LATTICES 



This part is devoted to the study of regular lattice folding in the absence of self- 
avoidance (phantom folding). Such foldings were considered for instance as discrete models 
for polymers (ID) or for "tethered" membranes (2D), equipped with an underlying regular 



elastic skeleton, such as crystalline or polymerized membranes ||2T| - |2^ . Our task is to first 
identify the folding degrees of freedom and to finally derive explicit statistical properties 
of large folded objects in the presence of bending rigidity. 



1. Phantom folding problems 

In all generality, we may define the phantom folding of a iD-dimensional lattice in 
a (i-dimensional target space with d > D a,s follows. Viewing the lattice as a simplicial 
complex with nodes, links, faces, ... I^-cells, a folding is a map from the lattice to IR'^ 
such that its restriction to each D-cell is an isometry, i.e. preserves the distance between 
any pair of points within the cell. The shapes of the D-cells are therefore preserved in 
a folding, but the relative position of two adjacent D-cells of the lattice may change as 
long as they remain adjacent to ensure a unique position of their common {D — l)-cell. In 
other words, the [D — l)-cell serves as a "hinge" between them. Such a folding is phantom 
in the sense that we allow cells to interpenetrate one-another, as well as distinct points 
of the lattice to occupy the same target position. A phantom folding is therefore defined 
by the corresponding folded configuration, without reference to an actual folding process 
which would continuously deform the original lattice. All the foldings considered in parts 
A-C below will be phantom. Moreover, we will restrict ourselves to the cases D = 1, d 
arbitrary or D = 2 and d = 2 or 3. We are mainly interested in the statistical properties 
of folded configurations of large portions of lattices for various ensembles corresponding 
to possibly restricted target spaces within IR'^. For instance, we may consider continuous 
foldings in which folding "angles" between adjacent D-cells vary continuously with specific 
distributions, or discrete foldings in which these angles only take finitely many values. For 
a good choice of these discrete angles, this simply corresponds to discretizing the target 
space into a d-dimensional lattice. 
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continuous 



discrete 



Fig.l: Folding of a segment of ID lattice in either continuous or discrete 
target space. The associated bending energy takes the form —KcosO for 
each pair of consecutive links with relative angle 9. An external force F may 
also be exerted at the ends. 

1.1. Folding: some trivial examples 

In this Section, we review two particularly simple folding problems, namely that of a 
ID lattice in any dimension, and that of the 2D square lattice in ci = 2. 

o ID folding 

As a preliminary exercise, let us study the case D = 1 of the phantom folding of a 
segment of length, say L of the lattice of relative integers Z. Viewing this segment of lattice 
as made of a chain of L + 1 nodes connected by L oriented links of unit length, a folding 
in IR'^ is entirely characterized, up to a global translation, by the sequence {ti}i<i<L 
of "tangent" vectors, all of unit length, which are the images of the links after folding, 
themselves connected into a possibly overlapping chain. We decide to attach to each 
folded configuration a weight 



w{{U}) = eM-E{{ti})) 



(1.1) 



where the energy of the configuration reads 



L-l 



L 




(1.2) 
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corresponding to a "bending" energy with bending rigidity K and an external force F 
between the extremities of the chain. We may consider several statistical ensembles for 
the tiS. Natural choices are (see Fig.l): 

(a) A continuum equidistributed set of unit vectors ti in IR'^. 

(b) A discrete equidistributed set with ti G {±ei, ±62, ■ • ■ , ie^} where {ei}i<i<d is the 
canonical basis of IR^. This choice alternatively amounts to requiring that the nodes 
of the chain have their folded images on that of a d-dimensional hypercubic lattice. 
Let us start with d = 1, in which case only the discrete ensemble (b) above makes 

sense. The model with tangent vectors ti = ±1 and weights (|1 . 1| ) ( p^ ) is nothing but the 
ID Ising model with spin coupling K and magnetic field F. The partition function on a 
segment of length L and with free boundary conditions reads: 



Zl{K,F)= {coshF+^^^^^^!^^^^^^ (e^cosh(F) + A)^-i 

/ , ^ e"^ + e'^sinh^(F)\ , ^ , 
+ coshF ^ (e^ cosh(F) - A)^"^ 



:i.3) 



with A = -^6-2^ +e2-^sinh^(F). We may define the order parameter Ai as the average 
orientation of the links through Ai = (l/L) ^^ti, which is nothing but the magnetization 
in the Ising language, with expectation value M = {Ai) = {l/L)dFl^og Zl{K, F). For 
large L and small F, we easily get 

M ~ e^^F (1.4) 

In particular, we recover the well known (ID Ising) result that there is no spontaneous 
orientation, i.e. M = for F for all finite values of K, expressing that the chain is 
always "crumpled" in the absence of external force. Still a "fiattening" transition occurs 
at K = +00 at which M = sign(F) = ±1 for F — 0, corresponding to a fully extended 
chain. 

This result is qualitatively the same for all target dimensions d. For instance, if the 
target space is the (i-dimensional hypercubic lattice (case (b) above), we get an average 
orientation At in the direction of the force F with expectation value: 



d-l + e^ 
d{d-l + e 



, , _^. F (1.5) 



for small F, which displays the same fiattening transition at K = +cx3. 
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Similarly, in the continuous case (a) above, we get 

'^ + 9iK) 

d{l-g{K)) ^ 

where g{K) = dK^og{j dHex.p{Kt • e)) for any fixed unit vector e. For instance, when 
d = 3, we have g{K) = coth(K) ~ 1/K. As g{K) tends to 1 only at K = +00, we find 
again a flattening transition at = +00. 



o Folding of the square lattice in d = 2 




Another particularly simple example is that of phantom folding of a rectangular piece, 
say of size Li x L2 of the 2D square lattice into a plane target. The image of each 
square cell is a unit square in the plane, entirely specified by the data of the two unit 
and orthogonal "tangent" vectors t^^^ and t^'^\ images of the south and west boundary 
links of the cell, pointing from their common origin (SW node). It is easily seen that 
preserving the shape and connectivity of the square cells implies that, up to a global 
rotation, t^^^ may take only the two values, say ±ei while t^^-* = ±62, and that moreover 
the value of t^^^ (rcsp. 

is the same for all south (horizontal) oriented links in a column 
(resp. all west - vertical - oriented links in a row). In other words, the folding constraints 
impose that the creases must propagate along straight lines all the way across the lattice. 
Each folded state is therefore entirely determined by the south (resp. west) ID tangent 
vectors along the southern (resp. western) border of the Li x L2 rectangle, namely by two 
independent sets of Ising spins {ti^^}i<i<Li and {tj^^}i<j<L2 (see Fig.2). In particular, the 
total number of folded states is 2^1+^2 ^ leading to a thermodynamic entropy per square 
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of s = linii;^ ^2^00 l/(-^i-^2)Log(2^i+^2^ = 0. The above factorization also allows to write 
the bending energy of the entire rectangle in terms of the spins as 

E, = -KL, - A-ii E if 42i (1.7) 

i=l j = l 

As in the ID case, we may also introduce a symmetry breaking field F coupled to the 
"orientation" of the folded state. For all 2D folding problems and for (i = 2 or 3, it is 
convenient to couple F to the projected area on a given plane of the target spacei. For 
(i = 2, this amounts to adding an energy 

= -^EE^S'^f ^ = -Fs,s2 (1.8) 

i=l 3=1 

where Sa = Ym=i ^I"^ total (algebraic) projected length in the direction a. The 

partition function may be written as 

Li L2 

Zl^^l,{K,F)= Y1 E e^'^^'^^ZL^{KL2,Si)ZL^{KL^,S2) (1.9) 

S]^ — L ]^ mod 2 S2—L/2 mod 2 

where Zl{K, S) is the coefficient of ^ in Zl{K, F) as given by Eq.(|1.3|), corresponding 
to the contribution of total projected length 5" to the ID partition function. In the ther- 
modynamic limit, it is convenient to introduce the free energy per square at fixed projected 
lengths Sa = (JaLa 

/(ai,a2) = lim --i-Log[e^^^^^'^i'^^ZL,(KL2, aiLi)ZL,(KLi, (T2L2)] (1.10) 
as well as the total free energy per square 

/= hm --i-Log[Zi„L,(K,F)]= min /(ai,a2) (1.11) 

Li,L2^oo 1ji1j2 0-1,0-2 6 [-1,1] 



For large Li and L2, the effective couplings KLi and KL2 in ( |1.9|) tend to infinity. For 
positive K, Zl^{KL2t Si = aiLi) is then dominated by two (symmetric) configurations 

^ Note that for a ID chain, the coupling to the "projected length" on a given line is equivalent 
to a force exerted at the ends of the chain in the direction of this line. 
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(a) (b) 

Fig.3: A typical configuration (a) minimizing or (b) maximizing the number 
of folds at fixed projected length Si. In case (a), the number of folds is at 
most one while in case (b), it is given by Li — IjSiI — Ss^fi- 

with at most one fold in the direction 1 (see Fig.3), hence, in the thermodynamic limit, 
we have 

f(ai,a2) = -2K-Faia2 (1.12) 

The minimum of this free energy is attained at ai = sign(F)(T2 = ±1 expressing the 
complete flatness of the lattice in the thermodynamic limit for any F. The corresponding 
total free energy reads 

f=-2K-\F\ (1.13) 

For negative K, Zl^{KL2, Si) is dominated by the configurations which maximize the 
number of folds in the direction 1 (see Fig.3), hence with exactly Li — — 5sj^,o such 
folds, contributing by an energy i^(l — 2|cri|) and by an entropy in the thermodynamic 
limit (the number of these configurations is clearly bounded by 2^^). This leads to a free 
energy 

/(cTi, a2) = 2K{1 -\ai\- jasl) - Fcnaa (1.14) 

The minimum of this free energy is now attained at cri = (72 = for |F| < —AK, and at 
CTi = sign(F)cr2 = ±1 otherwise, leading to 

f = m.m{-2K -\F\,2K) (1.15) 
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Fig.4: Phase diagram of the phantom folding of the square lattice in 2d, 
in the plane {K,F) of bending rigidity K and symmetry breaking field F. 
First order transition lines (thick lines) separate three regions with constant 
average projected area per square M = 0, +1, or —1. In each of these phases, 
the lattice is frozen in a unique state, either completely flat (M = ±1) or 
completely folded (M = 0). 

For K < and in the thermodynamic limit, the lattice is therefore totally flat for |F| > 
—AK, and folded for |F| < —AK. More precisely, the folded phase is dominated by a single 
conflguration with maximal number of folds in both directions, i.e. corresponds to the pure 
state of the completely folded lattice. The transition across the lines \F\ = —4K is flrst 
order, as well as that across F = for K > (see Fig.4). Note that, as opposed to the ID 
case, the flattening transition now occurs at = 0. 

Although it may seem at flrst rather pathological, the square lattice planar folding 
captures a number of essential features of regular lattice folding: 

- the non-local nature of the creases which propagate throughout the lattice; 

- the consequent absence of local excitations of the pure flat state; 

- the existence of a flrst order transition separating a completely flat phase made of this 
single flat state from a crumpled phase. 
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Rectangular Triangular 




Square-Diagonal Double-Triangular 

Fig. 5: The classification of two-dimensional compactly foldable lattices. 
1.2. Compactly foldable 2D lattices 

The details of folding are strongly lattice dependent. A good criterion for choosing 
the lattices to be folded is that they give rise to sufficiently many folded configurations. 
This is usually a prerequisite for allowing interesting (geometrical) phase transitions. For 
2D lattices, a basic requirement is to demand that the lattice might be completely folded 
onto one of its faces. Such a lattice will be called compactly foldable. This restricts very 
strongly the possible lattices: indeed, the corresponding lattice must have only one type 
of face, together with its finitely many possible rotations and refiections. This gives rise 
to a classification of the two-dimensional compactly foldable lattices: they fall into the 
four cases depicted in Fig. 5, namely rectangular, triangular, square- diagonal and double- 
triangular. 

The proof of this classification goes as follows . Let us concentrate on a node of the 
lattice. Each adjacent edge may serve as a hinge in the folding of its two adjacent faces, 
hence bisects the wedge made by its two neighboring edges around the node. Moreover, as 
the normal vectors alternate between neighboring faces in the completely folded state, the 
number of faces adjacent to every node is even, as is the number of edges adjacent to any 
node. Each node v is therefore the center of a regular star of say 2m„ edges forming angles 
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of -K/rriy, niy > 2, hence those angles are either right or acute. The faces are therefore 
polygons with at most four edges, and they can have four only if they are rectangles. This is 
the case where all m„ = 2, the rectangular lattice. Otherwise, all faces must be triangular, 
with right or acute angles. Such a face has angles 7r/mi, 7r/m2, vr/ma, with > 2, and 
^ l/m^ = 1. There are only three solutions up to permutation for (mi, m2, rris) namely 

(3, 3, 3) Triangular 

(2, 4, 4) ^ Square - Diagonal (1.16) 
(2, 3, 6) Double — Triangular 

This completes the proof. 

It is clear that the folding of the rectangular lattice is equivalent to that of the square 
lattice, in particular its 2d folding is trivial as discussed in the previous Section, with a 
vanishing thermodynamic entropy of folding. As it turns out, the three other cases share 
the property of having a non-zero thermodynamic entropy of folding per face 0,25|. In 



the following Sections, we will concentrate on the case of the triangular lattice, which we 
will fold both in (i = 2 and 3. 



2. Folding of the triangular lattice 
2.1. Generalities 

In this Section, we address the question of the folding of the triangular lattice. This 
is the simplest non trivial example of regular lattice folding and it captures all the generic 
properties of compactly foldable lattices. Another advantage of this model is that it bears 
links to some exactly solvable integrable model. Very generally, a folded configuration 
of the triangle lattice will consist in keeping the triangles equilateral of unit side, while 
the edges serve as hinges between neighboring triangles which may a priori form arbitrary 
angles. 

In the following, we will study two particular examples of folding: 

(a) The planar folding, i.e the case d = 2. 

(b) The folding on the Face Centered Cubic (FCC) lattice corresponding to d = 3 and 
suitably defined folding angles. 

The foldings of the triangular lattice may be characterized equivalently by link vari- 
ables, corresponding to tangent vectors, or face variables, corresponding to normal vectors. 
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Fig.6: Orientation of the links of the triangular lattice. The link (t^) and 
face (n/) variables. 

o Link variables 

We first start by orienting all the links of the triangular lattice in a coherent way 
so that all triangles pointing up are oriented counterclockwise and all triangles pointing 
down are oriented clockwise (see Fig.6). To each oriented link £, we associate its image 
in the folded configuration. To ensure that any elementary triangle has for its image an 
equilateral triangle of side 1 in JR!^, it is sufficient to require that all t^s are unit vectors, 
and that moreover 

te,+ti,+te,=0 (2.1) 
for the three oriented links £i, £2 and £3 around each triangle. 

o Face variables 

We can also describe a folding by face variables by considering the normal vector to 
any triangle /. For d < 3, it reads 

nf = e^h,Ate, (2.2) 

where ii and £2 are two consecutive oriented links around / and with e = +1 (resp. — 1) 
for triangles pointing up (resp. down) on the lattice. 

It is natural to think of this normal vector as a spin variable. The introduction of a 
bending energy amounts to a spin-like coupling K between normals of neighboring faces. 
As we shall see in the following, such energy is responsible for various conformational 
transitions. 

As for link variables, the face variables are not independent on each triangle. The six 
normal vectors for the six triangles around a node are correlated via a vertex-constraint. 
A more precise formulation of this constraint will be given in the case (a) of the planar 
folding, where it leads to an "11-vertex model", and in the case (b) of the folding on the 
FCC lattice, where it leads to a "96- vertex model" . 
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Fig. 7: The 11 possible folding environments of an elementary hexagon, cor- 
responding (from top to bottom) to : no fold, 1 fold (3 configurations), 2 
folds (6 configurations) or 3 folds. The pictorial representation on the left is 
only intended as a guide for the reader and only the final (phantom) folded 
state is relevant. The folds are indicated on the right by thick lines. 

2.2. Folding of the triangular lattice in d = 2: an 11 -vertex model 

The simplest case of folding, first introduced by Kantor and Jaric [^, corresponds 
to the planar folding, i.e. the folding in the plane M^. In this case, the angle between 
two neighboring triang lesi is either 0° (no fold), or 180° (complete fold). Folding an 
elementary hexagon leads to the 11 possible environments represented in Fig. 7. We recall 
that our definition of phantom folding is limited to the image of the lattice and does not 
distinguish between the different ways of reaching the same configuration. One can easily 
convince oneself that all folding constraints reduce to those of Fig. 7 around each elementary 
hexagon. This will become clear below in the language of tangent vectors. The problem of 
planar folding is therefore entirely characterized by the local rules of Fig. 7 on elementary 
hexagons. As such, planar folding is an 11-vertex model. 

We will take here as definition of tlie angle between two triangles tliat between their normal 
vectors. 
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Fig.8: The 11 possible spin environments (up to a global reversing of the 
spins) around a node. In all cases, the number of +1 (resp. —1) spins is a 
multiple of 3. 

o Face variables: constrained Ising spins 

The normal vector to each triangle may be described by an Ising spin which may take 
only two values a = +1 or a = —1 according to whether the triangle has been flipped 
or not in the folding. The folds are nothing but the domain walls for the spin variable. 
The 11 environments of Fig. 7 correspond to the 11 possible environments in terms of spins 
(up to a global reversal of the spins) indicated in Fig.8. A simple characterization of the 
allowed spin configurations is that the number of spins +1 must be equal to 0, 3 or 6, 
i.e. must be a multiple of 3, or equivalently: 

6 

J^a^ = mod 3 (2.3) 

for the six spins around each node. One easily checks that the number of such configura- 
tions is + (^) + (^) = 22 = 2 X 11 as it should. 

o Link-variables: the three-color model 

A natural question concerns the computation of the folding entropy which character- 
izes the exponential growth of the number Zjy of possible foldings as a function of the 
number of triangles (for a finite sub-lattice of the triangular lattice) . One defines the 
folding entropy s byll: 

s = lim —Log Zn = Log q . (2.4) 

A^— >+oo iV 

^ In this definition, one expects that the actual shape of the sub-lattice should not matter as 
long as the two dimensions scale in the same way and provided we take free boundary conditions. 
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Fig.9: Starting from the completely folded configuration, we may reverse the 
6 spins of some elementary hexagons (here in grey) and still get an allowed 
configuration provided the hexagons do not overlap (hard hexagons). 



One can interpret q as the number of degrees of freedom per triangle, with clearly, for 
constrained Ising spins, 1 < Q < 2. The existence of a non-vanishing entropy s (i.e. ij > 1) 
is not a priori obvious as, like in the square lattice case, there exists no configuration for 
which the folds are localized. A simple way of seeing this is to cut each configuration of 
Fig. 8 vertically in the middle and to notice that whenever a fold is present on the left half 
of the hexagon, then a fold is also present on the right side. The folds therefore propagate 
all the way from and to infinity from the left to the right. In other words, despite the 
locality of the rules of Fig. 7, the fiat configuration (without fold) has no local excitation. 
Numerically, Kantor and Jaric have shown that q ~ 1.21 and that there is therefore 
a non-vanishing folding entropy, as opposed to the square lattice case. This result may 
be understood by starting instead from the completely folded state, with a perfect anti- 
ferromagnetic order for the spins. We may then reverse globally all the spins of one single 
elementary hexagon and the configuration remains admissible (see Fig.9). This operation 
may be repeated for arbitrarily many hexagons as long as they do not overlap. There 
exist therefore local excitations of the completely folded configuration, leading to a non- 
vanishing entropy. 

The entropy may be computed exactly in the link- variable picture [Q. Starting from 
a given triangle, the vectors ti^, i = 1, 2, 3 around this triangle are three unit vectors in IR^ 
with vanishing sum, therefore of the form ti-^ = A, ti^ = B = IZA and ti^ = C = IZ^A = 
TZ~^A where A is a unit vector and TZ the rotation of angle 120° (or —120°). One can then 
easily convince oneself that the vectors ti for all the links of the lattice may only take these 
same three values A, B or C, with moreover the constraint that the three "colors" A, B 
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Fig. 10: The three allowed values (colors) for ti are A, B = TZA and C = 
TZ~^A. The relation between link variables ti and the face variables ai = ±1 
around a node reads ti^i = TZ^Hi. The link variable is well defined if = ti, 
that is if fOI) is satisfied. 



and C must be present once and only once around each triangle [Q. A planar folding of 
the triangular lattice is therefore equivalent to a tricoloring of the links of the lattice with 
three colors which must be distinct around each triangle. The color simply represents the 
orientation of the link in the folded configuration, among three allowed values. 

In this language, the spin of a face tells us about the cyclic order of the colors around 
this face with a = +1 (resp. cr = —1) if the colors are in the order A,B,C (resp. A,C,B) 
in counterclockwise direction. In particular (see Fig. 10), the variables t^, z = 1, ■ ■ ■ , 6 of 
the links around a node are related to the spins by ti^i = TZ^Hi. As 7^3 = 1, the link 
variable is well defined after one turn around the node if and only if the relation ( p.3| ) is 
satisfied. In this language, all folding constraints ensuring that the link variable is well 
defined along any closed loop clearly reduce to Eq.( |2.3|) around each elementary hexagon. 

In the dual language, a planar folding therefore corresponds to the coloring of the 
links of the hexagonal lattice with three colors, such that no two adjacent links be of the 



same color. The entropy of this three-color problem was obtained by Baxter [^, and 
can therefore be reinterpreted as the planar folding entropy per triangle of the triangular 
lattice, with the result: 



■i, v/3n(3n - 2) r(2/3) 2ti ^ ' > 



n 



(2.5) 



Its numerical value q = 1.208717... (s = Logg = 0.189560...) is in very good agreement 



with the estimate of Ref. [E6 . 



18 



F 







u.o 












0.4 




M=+l 








.0.2 




(gap 
















-0.3 


-0.1 

M=0 


> 


/o.i 


0.3 




Cjolded^ / 


-0.2 




M=-l 








-0.4 




(tiat) 








-0.6 









Fig. 11: Phase diagram for tlie planar folding in the presence of a bending 
energy term Ei, = —Kcos{9) and a symmetry breaking field F. A first 
order transition at Kc ~ 0.1 separates the crumpled phase where the average 
projected area M = (a) = from a fiat phase where M = 1 (or —1). This 
latter phase is frozen in the pure state of the completely fiat surface. Within 
the crumpled phase, a continuous piling-up transition now occurs at Kp ~ 
—0.28 []^,^^ between a disordered phase with M^t = and an ordered phase 
with Mst ^ 0. 

We may finally consider the phase diagram obtained in the presence of a bending 
energy Ef, = —Kcos{9) per link with folding angle ^ = 0° (no fold) or 180° (complete fold). 
Like in the square lattice case, we also introduce a symmetry breaking field F conjugate 
to the projected area {Ep = —Fa per face) and which plays the role of a lateral tension 
applied to the surface. The resulting phase diagram, depicted in Fig. 11, was obtained by 
a transfer matrix approach ^ and later corroborated by a variational study [p8| . 

Comparing this result to that of the square lattice (Fig. 4), we note a qualitative 
resemblance, with three first order transition lines separating regions with M = 0, ±1. 
There are however two important differences. Firstly, the phase M = is no longer frozen 
in a single state due to the existence of elementary excitations of the fully folded state. 
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More precisely, let us introduce the "piling- up" (staggered) order parameter Mgt = {{M.a " 
M.s])/2) where M./^ (resp. M.s]) is the average projected area of the triangles pointing up 
(resp. down) in the original triangular lattice. A continuous piling-up transition now occurs 
at a negative value ~ —0.28 above which the crumpled phase is disordered (Mgt = 0) 
and below which the lattice tends to be piled- up onto a single triangle (Mgt 7^ 0). Note 
that the transition lines between the crumpled and flat phases are no longer straight lines. 
Secondly, the triple point where the three transition lines meet is at a positive value Kc ~ 
0.1. This shows the existence of a first order flattening or crumpling transition between a 
crumpled phase for K < Kc and a completely flat phase for K > Kc corresponding again 
to a pure state without any fold, with M = 1 (or —1). The existence of such a crumpling 
transition was also predicted in continuous models of tethered membranes [pTj - plf . The 
first order nature of the present transition seems however to be due to the discrete (spin- 
like) nature of the (normal vector) degrees of freedom. For K = 0, comparing the energy 
per triangle —F in the fiat state (no entropy) and the entropy Logq in the folded state 
(no energy since M = 0), we deduce that the transition lines cross the vertical axis at 
values F = ±Logq ~ ±0.189. The existence of a strictly positive value of the crumpling 
transition rigidity is therefore directly linked to the existence of a non-zero entropy of 
folding. In some sense, the critical point i^T = of the square lattice folding has been split 
into two transition points Kp < < Kc (details may be found in Refs. [[7|,|8|,p8| ) . 

2.3. Discrete folding of the triangular lattice in d = 3.- a 96-vertex model 

A second simple example of folding of the triangular lattice consists in picking the Face 
Centered Cubic (FCC) lattice for target space 0. The latter may be seen as a discretization 
of IR"^ for which the faces of all the elementary 2-cells are equilateral triangles. More 
precisely, the FCC lattice may be viewed as a regular piling-up of elementary octahedra, 
completed by tetrahedra, as indicated in Fig. 12. 

The folding of the triangular lattice on the FCC lattice authorizes four folding angles: 
0° (no fold) 180° (complete fold) 109°28' (acute fold, i.e. on the same tetrahedron) and 
70°32' (obtuse fold, i.e on the same octahedron). The link variables ti then take their 
values in a set of 12 allowed values (ie^ ±6-,)/ \/2 with 1 < ^ < j < 3. These 12 vectors are 
in one-to-one correspondence with the oriented edges of an elementary octahedron of side 
1 (see Fig. 13). A simple way to label them is to first color them in three colors A,B,C such 
that all four vectors in a given plane of the octahedron of Fig. 13 be of the same color. We 
then complete the color by a complementary variable defined as follows: for each vector of 
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Fig. 12: The FCC lattice viewed as a piling-up of octahedra completed by 
tetrahedra. 

color A or B, we define a variable C = ±1 indicating on which side of the plane of color C 
on the elementary octahedron this vector lies (see Fig. 13). Similarly, we define a variable 
S = ±1 for the vectors colored A or C and a variable ^ = ±1 for those colored B or 
C. Each of the 12 unit vectors is entirely specified by its color and the value of the two 
complementary variables. 



Fig. 14: The 6 link variables and the 6 pairs of spins around a node. 

The three vectors t^. , z = 1, 2, 3 around a given face must have a vanishing sum: they 
must be chosen among the 8 triplets of unit vectors with vanishing sum (in correspondence 
with the 8 faces of the octahedron of Fig. 13), which, with the 3! permutations of the three 
vectors, leads to 48 possibilities. For each triplet, the three vectors have different colors, 
and have the same values of A, B and C, i.e are of the form {A, B, •), {A, •, C), (•, B, C). 
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Fig.13: The 12 oriented edges of an elementary octahedron of side 1 provide 
a set of 12 unit vectors for the aUowed hnk variables of the FCC folding of 
the triangular lattice. We attach colors A,B,C to these vectors as indicated. 
We have three color planes, each containing four vectors. We further attach 
variables A, B and C G {±1}- Each vector is coded by a triplet of the form 
{Ai B, ■), {A, ■,€) or B,C), where the position of the missing variable (coded 
by a dot) corresponds to the color of the vector at hand. Note that the three 
vectors around each face share the same A, B, C variables. 

We may then follow the 6 link variables ti around a node in clockwise direction (see 
Fig. 14). Starting from ti which takes one of the 12 allowed values, the vector t2 may take 
only 48/12 = 4 possible values coded by two face spin variables cr = ±1 and z = ±1. The 
variable a indicates the change of color from ti to t2, with cr = +1 if the color "increases" 
in cyclic order ABC and a — —1 otherwise. The variable z is simply the product ABC 
corresponding to the face {ti,t2, —ti — ^2} on the octahedron of Fig.13. 

Going back to the original folding problem, we may now characterize entirely the 
relative folding state of any two neighboring triangles by the relative values of (ui, zi) vs 
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[(^2tZ2) on these two faces. We have the foUowing correspondence 



Z2/ZI 


0-2/0-1 


type of fold 


1 


1 


no fold 


1 


-1 


complete fold 


-1 


1 


acute fold 


-1 


-1 


obtuse fold 



After a complete turn around a node, the constraint that we recover the same tangent 
vector translates into the following two constraints on the spins. The first constraint 



o"i = mod 3 



(2.6) 



is the same as the constraint ( |2.3[ ) for planar folding, and ensures that we recover the same 
color of the tangent vector. The second constraint deals with the z variableJi and reads 



n 



i such that: 

X mod 3 



X 



1,2,3 



(2.7) 



The left hand side simply counts the number of sign changes of the z variable across the 
links of color A, B or C (one color for each value of x), which are nothing but the sign 
changes of the complementary variable A, B and C respectively. Eq. ( p.7|) ensures that 
the we recover the same value of the complementary variables after a complete turn. 

In terms of constrained spins, it is easy to count all the environments allowed by the 
constraints (|2.6|) and ( p.7|) around a vertex. There are 96 such environments, which are 
represented in Fig. 15. As such, the FCC folding of the triangular lattice is a 96- vertex 
model. 

It is interesting to note that this model has the following three sub-models (see Fig. 15): 
1. An 11-vertex model corresponding to the planar folding obtained by keeping only 
either unfolded or completely folded links. This amount to requiring that z = 1 for 
all the triangles, in which case is automatically satisfied and one recovers the 

constraint 



Note that the relation for a; = 3 is a consequence of those at x = 1,2. 
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plane 
(11 vertices) 



octahedron 
(16 vertices) 




tetrahedron 
(11 vertices) 



96 vertices 



Fig.15: The 96 folding configurations of an elementary hexagon. The com- 
plete folds are indicated by thick lines, the obtuse folds by thin lines and 
the acute folds by dashed lines. For each environment, we have indicated its 
multiplicity corresponding to rotations of the hexagon. We also indicated the 
vertices which have to be retained for the three sub-models corresponding to 
the planar folding ( 1 1 vertex) , the folding on a single tetrahedron ( 1 1 vertex) 
and that on a single octahedron (16 vertex). 

2. An 11-vertex corresponding to the folding on a single target tetrahedron, obtained by 
keeping only acute or complete folds. This amounts to requiring that z = ea where 
e is 1 on triangles pointing up and —1 on triangles pointing down. Once again, the 
relation (|2.7|) is automatically satisfied. 

3. A 16-vertex model corresponding to the folding on a single target octahedron, obtained 
by keeping only obtuse or complete folds. This amounts to having a perfect anti- 
ferromagnetic order a = e, with e as defined above, in which case the relation ( |2.7| ) 

becomes HLi ^i^i+i = 1- 
Note that the constraints ( p.6| ) and (|2.71 ) are invariant under the change (cr, z) (a, eza). 
One thus gets a duality relation which consists in exchanging globally all the unfolded 
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edges with those carrying an acute fold, as apparent in Fig. 15. The sub-model 3 above 
is self-dual while the two sub-models 1 and 2 are dual of one another. In particular, the 
entropy of folding on a single tetrahedron is the same as that for the planar folding, and 
given by Eq. (|2.5|) . 

The entropy of folding on the FCC lattice is not know so far. Numerical estimates 
show that q ~ 1.43(1). We may easily show that q > \/2 = 1.414- ■ ■ by estimating 
the entropy of the sub-model 3 of folding on a single octahedron. Indeed, the constraint 
11^=1 ^i^i+i = 1 amounts to requiring that z is the product rjy-^^rjy^rjy^ on the three nodes Vi, 
f 2 and Vs adjacent to the face at hand of a node variable ry^ equal to ±1 and independent 
on each node of the lattice. The model 3 has therefore an entropy of 2 per node, i.e. \/2 
per triangle. 
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-1.0 -0.8 -0.6 -0.4 -0.2 0.0 

K 

Fig. 16: The proportion of the different types of folds as a function of the 
rigidity modulus K (energy E{, — —Kcos{9)) in the range K < 0. 

We may finally consider the phase diagram obtained by introducing a bending energy, 
again in the form of a link energy Ef, = —Kcos{9) where 9 is the folding angle (none. 
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complete, acute or obtuse). This phase diagram was obtained in Ref.0 by use of a vari- 
ational method (see also Ref. [^). As for planar folding, one finds in the regime K > a, 
first order fiattening/crumpling transition between a crumpled phase for small values of K 
and a completely flat phase for large values of K. The situation for K < (i.e. for which 
folds are favored) is richer: the figure 16 displays the proportion of the different types of 
folds as a function of K. One clearly distinguishes two successive folding transitions. A 
first, discontinuous transition at -ftT ~ —0.3 separates a regime where the folding occurs 
preferentially on octahedra from a regime where it occurs preferentially on tetrahedra. A 
second, continuous transition at -ftT ~ —.85 separates this last regime from a regime where 
the folding is essentially maximal with dominance of complete folds. 
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Fig. 17: The three order parameters O, T and P as defined in the text as 
a function of the rigidity modulus K (energy Eb — —Kcos{9)) in the range 
K <0. 
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These transitions may also be viewed on the three order parameters O, T and P 
defined as 

T={zea) (2.8) 
P^{z) 

and represented in Fig. 17. The interpretation of these order parameters is that a non- 
vanishing value of O (respectively of T) shows a tendency of the lattice to fold preferentially 
around octahedra (respectively around tetrahedra). Similarly, a non- vanishing value of P 
shows a tendency of the lattice to preferentially remain within planes of the FCC lattice. 

This completes our study of phantom folding of 2D regular lattices. Here we have 
concentrated on the folding of the triangular lattice in d = 2 and 3 dimensions. We expect 
the same type of conformational transitions to take place for higher target dimensions d, 
with a first order crumpling transition to a completely fiat state at a positive value of the 
bending rigidity K, a continuous (Ising type) piling-up transition at a negative value of K 
toward a limiting completely folded state, and a fine structure of intermediate transitions 
of wrapping around intermediate-dimensional regular solids. Finally, we expect a similar 
scenario to happen for all the compactly foldable lattices of Fig. 5 (including the rectangular 
lattice degenerate limit with zero entropy) . This scenario is confirmed in the case 

of the square-diagonal lattice in Ref . [p5 1 . 
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PART B: LOOP MODELS ON REGULAR LATTICES 



In this part, we show that folding problems are particular instances of a larger class 
of problems, so-called fully packed loop models with possible reformulations as height 
models for which the loops play the role of contour plots. We focus here again on regular 
lattices and illustrate the techniques of solution in the case of fully packed loops on the 
hexagonal lattice. These include Bethe Ansatz calculations and effective Coulomb gas 
descriptions |p2 . 



3. Loop gas and height model reformulations of the triangular lattice folding 

Very generally, all lattice folding problems can be reformulated as fully packed loop 
gases on the dual lattice. We will illustrate this property in the two cases of the d = 2 
(planar) and d = 3 (FCC) foldings described in the previous Section. 

3.1. Fully packed loop gas formulations of the 2d folding 




Fig. 18: Example of tricoloring of the links of the hexagonal lattice and the 
associated configuration of oriented fully packed loops. 

Let us first consider the planar folding of the triangular lattice. As seen in Section 2.2, 
it is equivalent to the tricoloring of the links of the hexagonal lattice with three colors A, B 
and C, required to be distinct on the three links adjacent to any node of the lattice. Note 
that the hexagonal lattice is bipartite, i.e. it may be decomposed into two sub-lattices of 
black, resp. white nodes, such that the neighbors of a black node are white and conversely. 
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The black or white color of a node is equivalent to the "pointing up" or "pointing down" 
nature of the dual triangle. Ignoring the links of color A, the links of color B or C form 
loops on the hexagonal lattice along which the colors B and C alternate (see Fig. 18). By 
orienting all the links of color B from the adjacent black node to the adjacent white node, 
and the links of color C from the adjacent white node to the adjacent black node, we 
orient each loop in a coherent way. Changing the orientation of a loop simply amounts to 
interchanging the colors B and C along this loop. The gas of oriented loops thus obtained 
is self-avoiding and fully packed, in the sense that each node of the lattice is visited by a 
loop. 

Conversely, given any set of fully packed oriented loops, it is possible reproduce the 
associated colors for all the links, and therefore rebuild the associated folding configuration. 
This shows that the planar folding of the triangular lattice is equivalent to a gas of fully 
packed oriented loops on the hexagonal lattice. 

The orientation of the loops may be rephrased into a weight n = 2 per unoriented 
loop. More generally, we may consider the so-called FPL{n) model of a gas of fully packed, 
self-avoiding and unoriented loops on the hexagonal lattice with a weight n per loop pO[] . 
This is to be compared with the loop gas formulation of the so-called 0{n) model on the 
hexagonal lattice which also amounts to consider self-avoiding loops with a weight n per 
loop but without the requirement that every node of the lattice be visited by a loop ||3^ . 

Note that the 0{n) model is originally defined as a spin model with 0{n) symmetry. 
Similarly, the FPL{2) model alternatively describes the groundstates of the antiferromag- 



netic XY model (spins with 0{2) symmetry) on the Kagome lattice [^^. Indeed, the nodes 
of the Kagome lattice are in one-to-one correspondence with the edges of the triangular 
lattice. The tangent vectors of the folding configuration yield spins on the nodes of the 
Kagome lattice. The constraint ( p.l|) of vanishing sum around a triangle is nothing but the 
condition for the minimization of the antiferromagnetic energy t^-^ - tg^ +^^2 '^^s +^^3 '^^i (to- 
gether with \tg. \ = 1) around each elementary triangle of the Kagome lattice. In the same 
spirit, the FPL{1) model describes the antiferromagnetic groundstates of the Ising model 
on the triangular lattice, as illustrated in Fig.l9-(a,b). The loops on the (dual) hexagonal 
lattice are easily identified with the links dual to the -f|— (i.e. energetically favored) links 
on the triangular lattice. The fully packed requirement amounts to maximizing the number 
of these favored links (exactly 2 per triangle). 
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Fig. 19: Equivalence between (a) antiferromagnetic groundstates of tiie Ising 
model on the triangular lattice, (b) FPL{1) loop configurations on the dual 
hexagonal lattice, and (c) Solid-On-Solid (SOS) interfaces for piling-ups of 
cubes viewed in perspective. In (a) , there is exactly one frustrated link (with 
equal adjacent spins) per triangle (dashed lines). The frustrated links are dual 
to the unoccupied edges in (b) while the links dual to the -(-|— links form a 
gas of fully packed loops on the hexagonal lattice. Erasing the frustrated 
links yields a rhombus tiling of the plane, also interpreted as a piling-up of 
cubes viewed in perspective (c). 

3.2. Fully packed loop gas formulations of the 3d folding 

In the case of the folding of the triangular lattice on the FCC lattice, we have seen 
that the link variables again induce a tricoloring of the links of the triangular lattice, 
and thus, by duality, of the links of the hexagonal lattice. This coloring by three colors 
A, B and C, as we just saw, is equivalent to the gas of the fully packed BC loops. We 
may consider as well the two other systems of fully packed loops made of the links AC or 
the links AB. The three systems of loops are strongly correlated as the knowledge of one 
of them allows us to reconstruct the two others. The colors are not sufficient however to 
characterize entirely the original 3d folding. To this end, we have introduced the additional 
link variables A,B,C = ±1 indicating, in the elementary octahedron, the side of the plane 
of color A,B,C respectively on which the link variable at hand lies (see Fig. 13). If we now 
consider any BC loop, it is easy to see that the value of A is constant along the loop. 
Indeed, a change of sign of A requires to cross the plane of color A on the elementary 
octahedron and therefore to pass through a link of color A, which in turn requires to pass 



30 



from one BC loop to another one. Moreover, the value of A is independent on each of 
the BC loops. The same is true for the variable B on the AC loops and the variable C on 
the AB loops. As a consequence, we may represent any folding on the FCC lattice by a 
tricoloring of the hexagonal lattice, completed by a ^2 spin variable on each of the BC, 
AC and AB loops. The enumeration of 3d foldings is therefore performed by counting 
tricolored configurations with a weight 2 for each of the AB, BC and AC loops. 

To conclude, we may write the partition functions for 2d and 3d foldings as partition 
functions for edge-tricolorings or fully packed loop g 

^plane= ^ = E 2# ^^^^^ 

Tricolorings Fully packed loops 

Zpcc = E 2* ^°°P' X 2* ^""P"" X 2* ^°°P' 

Tricolorings 



Other folding problems were considered , such as the planar folding of the square- 
diagonal lattice, or of the triangular-diagonal lattice. All these problems are equivalent to 
gases of fully packed loops on appropriate lattices. 

3.3. Height models 

We shall now transform the FPL(2) configurations into those of a height model defined 
as follows P^. We start with a gas of oriented fully packed self-avoiding loops. To each 



face of the hexagonal lattice, we associate a "height" X in M'^ whose variation from one 
face to another depends on whether the separating edge is occupied or not by a loop. 
The precise rules for the variation of heights are indicated in Fig. 20 and use explicitly the 
bipartite nature of the hexagonal lattice, with different conventions according to whether 
the link is oriented from a white or from a black node. These rules use three elementary 
height differences A, B and C. If we now replace the oriented loops by the corresponding 
tricoloring of the links as in Fig. 18, the height difference is simply given by the color A,B 
or C of the crossed link. 

Inspecting the allowed configurations around a node, we see that we must impose the 
constraint A -\- B -\- C = Q m. order to have a well-defined height after one turn around a 
node (see Fig.20-(a)). This naturally leads to a two-dimensionalYvei^i X. As we shall see 
below, the components of X become two scalar fields in a continuous effective description, 
and the FPL (2) model is indeed a critical model described in the continuum by a conformal 
field theory (CFT) of central charge 2 corresponding to these two height components ||35|| . 
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Fig. 20: Definition of the height variables X for fully packed loops on the 
hexagonal lattice. The rule explicitly uses the bicoloring of the nodes of the 
hexagonal lattice. In the case of fully packed loops, the only possible environ- 
ment is (a) (up to obvious symmetries), and the consistency of the definition 
of the heights imposes A + 5 + C = 0. In the case of loops not necessarily 
fully packed, the presence of unvisited nodes with the environment (b) im- 
poses moreover that ^ = 0, in which case B = —C (= 1 for instance) and 
the rules become insensitive to the bicoloring of the nodes of the hexagonal 
lattice. 

A symmetric solution consists in choosing A to be a unit vector, B = TZA and C = TZ^A 
where TZ is the rotation of angle 120° in the plane. With this choice, for any given oriented 
loop configuration, we may identify the height X of a face with the position of its dual 
node in the associated folded configuration of the triangular lattice in the plane. The above 
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symmetric solution will be used from now on. 

As described above, the 0{2) model corresponds to releasing the full-packing con- 
straint by allowing for nodes which are not visited by loops. In this case, we must also 
include the configuration of Fig. 20- (b) corresponding to an unvisited node. The consis- 
tency of heights after one turn around the node now imposes the extra constraint A = 0, 
and thus B + C = 0. This now leads to a one- dimensional height. The corresponding 
loop gas is critical for a suitable choice of edge weight. At this critical point ||3^, the 
model is described in the continuum by a CFT of central charge 1. Notice that the rules 
defining the heights in this case are now insensitive to the bipartite nature of the lattice 
(see Fig.20-(b)). 

Surprisingly, the full-packing constraint infiuences the value of the central charge, and 
thus the universality class of the model. More precisely, as first recognized in Ref. , the 
net effect of the full-packing constraint has been to increase the central charge by 1. We 
shall recover a similar property below in the case of random lattices. 

For n = 1, the FPL{1) model may be described by a one- dimensional height. This is 
readily seen on Fig.l9-(c). Indeed, any configuration of the model may be interpreted as 
a piling of unit cubes in IR"^ viewed in perspective from the direction (1,1,1) and whose 
free surface is connected and has no overhangs. This surface may be entirely described 
by a single height h = X ■ A with X obtained from the rules of Fig. 20 for some arbitrary 
orientation of the loops (indeed this height is insensitive to the choice of orientation as 
B ■ A = C ■ A). As we shall see in the next Section, the model is critical with central charge 
1, as opposed the case of the 0(1) model which has central charge (in its dense phase). 

4. Exact solutions via Bethe Ansatz: the example of the FPL{1) model 

Exact solutions for the general FPL{n) models with arbitrary n may be obtained by 
transfer matrix methods together with Bethe Ansatz type assumptions on the eigenvectors. 
The existence of such solutions is granted by the integrability of the models. The Bethe 
Ansatz needed in the general solution is of "nested" type, which makes the explicit solutions 



quite technical [^^. In the case n = 1 however, the problem simplifies drastically and 
reduces to a free fermion model. The latter may be treated via a simple coordinate Bethe 
Ansatz. For pedagogical purposes, we choose to describe only this case in detail in the 
following Sections. In particular, we will compute the value c = 1 of the central charge 
for the conformal field theory that describes the large distance behavior of the model. We 
will return to the general case in Section 5 where we present an effective field theoretical 
description of FPL{n) via Coulomb gas. 
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4.1. FPL{1) and rhombus tiling of the plane 

In this Section, we address the FPL{1) model, which, as iUustrated in Fig. 19, may be 
rephrased as the problem of tiling a domain of the plane by means of any number of the 
three following rhombic tiles: 



a = / / P = \ \ 7 = < ) (4.1 




with edges of unit length and angles of 60° and 120°, all weighted by 1. As it turns out, the 
specific shape of the domain to be tiled is extremely important. One can show for instance 
that a large rhombus of size N x M may be tiled in exactly (^^^) distinct manners, 
resulting in a vanishing thermodynamic tiling entropy s = lim^r m^oo Wm^'^s{^~n^) = 0- 
On the other hand, a large hexagon with size N x M x P may be tiled in many more ways, 
namely ^6 



N M P 

i=ij=i k=i 

which results in a non-vanishing thermodynamic entropy of tiling. It is not our purpose 
here to study the influence of boundary conditions on the entropy of tiling, but just to 
compute it in some generic situation, namely that of a large cylinder represented by a 
rhombus of width L and length M ^ oo, the two longitudinal sides of the rhombus being 
glued. 

In the next Sections, we use the following strategy, quite standard for solving integrable 
lattice models (the particular example at hand is borrowed from a course delivered by B. 
Nienhuis at the 1997 Beg-Rohu school "Ecole de physique de la matiere condensee"). 
We flrst deflne the row-to-row transfer matrix Tl of the model, with periodic boundary 
conditions along a row of L tiles. The partition function of the model on a rhombus of size 
L X M , say with doubly periodic boundary conditions (torus) reads then 

Zlm = IV(Tf ) (4.3) 

To access the thermodynamic properties of the model, we must take a large L, M limit. 
Concentrating on a strip of flxed width L, the thermodynamics (M oo) is simply 
governed by the largest eigenvalue of T^, with a thermodynamic tiling entropy per tile of 

s= lim --\-Log Zl,m (4.4) 

We then use the expected conformal limit of the model to identify the central charge by 
studying the flnite size effects on a cylinder of flnite but large width L. 
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4-2. Transfer matrix and eigenvalue equations 



Consider a horizontal row, tiled with a number of tiles a, (3^ 7, and with total length 
L. For simplicity, we deform slightly the three tiles into the following 



a 



P 



1 



(4.5) 



This amounts to consider the partition function of the rhombus tiling model on a rectangle 
of size LxM with its vertical edges identified so as to form a cylinder. The transfer between 
a tiled horizontal row to the next (say from bottom to top) is then described by assigning 
an occupation number = 0, 1 (resp. n[) to each lower (resp. upper) horizontal edge, 
resulting in the transfer matrix element (TL){n'^,n^,...,n^},{ni,n2,...,ni,} from a configuration 
n of occupation numbers of a horizontal row of edges to the next n', that typically reads 
pictorially 



(rL){l,0,l,0,... ,0,1, 1},{1,1,0,0,..., 1,0,1} 

10 10 







• • • • 











(4.6) 

Note that the number /c of I's is conserved from a row to the next, and counts the total 
number of a or tiles in the row. 

Starting from a row with only 7 tiles, the next row can only consist of 7 tiles as well. 
This forms the "fundamental" invariant sector, with A; = 0. This unique possibility is 
summarized by the transfer matrix element 

(rL){o,0...,0},{0,0...,0} = 

(4.7) 

Note that all the upper and lower horizontal edges are unoccupied. We conclude that the 
"no-occupied" edge state is an eigenstate of Tl with eigenvalue Aq = 1. 

An excitation of this groundstate is obtained by replacing one pair of adjacent halves 
of 7 tiles along the row by either an a or a /? tile. Assume the occupied lower edge of the 
tile is in position x & Zl (with periodic boundary conditions L = 0) along the row, then 
the upper one is necessarily in position x (a tile) or a; + 1 {P tile). Let us denote by v{x) 
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the lower edge configuration vector {0, 0, 0, 1, 0, 0} with a 1 in position x. We then 
have 

Tlv{x)= v(x) + v{x+l) 

X x+1 




(4.8) 



X X 



We must consider multiple excitations of the groundstate as well. In the /c-excitation 
sector, k of the pairs of adjacent 7 half-tiles must be replaced by a or tiles. This results 
in the appearance of k occupied lower and upper edges. This conservation of the number 
of occupied edges from one row to the next is an important feature of the model. It allows 
to break the action of the transfer matrix Tl into L + 1 sectors, according to the number 
of excitations k e {0, 1, 2..., L}. The invariant space corresponding to k excitations has 
dimension (^) . For instance, in the case of k = 2 excitations, we must replace two of the 
pairs of adjacent 7 half-tiles by a's or /3's. The corresponding lower occupied edges are in 
positions say x and y,0<x<y<L — 1. Let us denote by ?;(,x, y) the edge configuration 
vector with only O's except for two I's in positions x and y. Just like in the 1-excitation 
case, the occupied upper edges must be in positions x or x + 1 and y or y+1, with moreover 
the constraint that they be distinct. Forgetting about this constraint, we would simply 
have 

Tl v{x, y) = v{x, y) + v{x + 1, y) + v{x, + 1) + v{x + 1, |/ + 1) (4.9) 

Imposing the constraint on the right hand side requires that, whenever y = x + 1, the 
second term should be removed, as well as the third one whenever y = L — 1 and x = 0. 
There is a much simpler way of automatically ensuring this: we simply enhance the domain 
of definition to a; e [0, L] and demand that v{x, x) ~ for all x, and that v{x, L) = v{0, x) 
for all x. 

The fc-excitation sector is easily constructed. Let v{xi, X2, Xk) be the edge- 
configuration vector with I's in positions xi < X2 <■■■<: Xk, with the convention that it 
vanishes whenever any two consecutive positions are equal, and that v{xi,X2, Xk-i, L) = 
v(0,a;i, ...,a;fc_i). Let us also introduce the translation operator Tj that adds 1 to the j-th. 
argument of such a vector, namely 

Tjv{xi, X2, a;j_i, Xj,Xj+i, Xk) = v{xi, X2, Xj^i, Xj + 1, Xj+i, Xk) (4.10) 



36 



Note that the result is the zero vector when Xj^i = Xj or xj + 1. The transfer matrix then 
acts in the /c-excitation sector as 

k 

Tlv{xi, ...,Xk) = + rj) v{xi,...,Xk) (4.11) 

where / is the identity and the various translation operators commute with one-another. 
The vanishing constraint on the vectors takes care of all the steric constraints for the a 
and /3 tiles. 

We now look for eigenvectors of in the form 

Yl f{xu---.Xk)v{xu---,Xk) (4.12) 

0<xi<X2<---a;fc<L — 1 

in the /c-excitation sector. The corresponding eigenvalue equation reads 

k 

A^yC^) = J](/ + r,) y^'^) (4.13) 

3=1 

This equation must be solved within the invariant -dimensional vector space of the /c-th 
sector, with the boundary conditions 

/(..., ...) = 

(4.14) 

/(xi,a;2, ...,Xfc-i,L) = /(O, xi, Xfc-i) 

4-3. Bethe Ansatz 

In view of the equations ( [4.13| ), we may naturally think of trying some simple Ansatz 
for the candidate eigenvectors, based on the eigenvectors of the translation operator r. In- 
deed, for /c = 1, f{x) = produces an eigenvector for r, for any non-zero z. The /c- variable 
version is /(xi,...,Xfc) = z^^ Z2'^ ....z^'^ , for some non-zero complex numbers zi, Z2, Zk- 
This Ansatz would work perfectly except that it violates the boundary condition, that / 
vanishes whenever two consecutive arguments are equal. 

The next fundamental remark is that if we still ignore the boundary condition, any per- 
mutation a G S'fc of the x's would yield an equally acceptable candidate /cr(xi, X2, Xk) = 
^x^(i)^x^(2) ^^^^^x^(fe)^ with the same eigenvalue. 

The idea of the Bethe Ansatz is then to combine these two ideas and look for a solution 



of the eigenvalue equation (|4.13| ) that is a linear combination of the /cr, a G S'fc, and that 
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incorporates the boundary conditions ( [4.14|) . In the present problem, the answer is unique 
and reads 

f{xi,X2,...,Xk) = det{z^')i<ij<k (4.15) 



Moreover, the periodic boundary conditions (second line of ( [1.14| )) along the strip are 
satisfied iff 

= (-1)'=-^ for ah i = 1,2,..., k 

(4.16) 

Zr ^ Zs for all r, s = 1,2, k, r^s 

where the extra condition that the Zi must be distinct ensures that the eigenvectors are 
non-zero. These are the celebrated Bethe Ansatz equations, in our particularly simple 
case. Note that we indeed get a basis of the -dimensional invariant space by taking / 
as in Eq. (^4.15|) with 

''-J 



= g2^^(n,-V)/^' j = 1^2,..., k 

' (4.17) 
< ni < n2 < ... < Uk < L — 1 



The corresponding eigenvalues read 



k 

k~l 



Akin^, ...,nk) = Hil + e^-K"— )/^) (4.18) 

3 = 1 

4.4- Continuum limit: largest eigenvalue 

The continuum thermodynamic limit of the model corresponds to L, M oo. For L 
finite and M ^ oo, the thermodynamic limit of the free energy of the model is governed 
by the largest eigenvalues of Tl. More precisely, we may compute the entropy of tiling 
per row of a strip of infinite length and finite width L by taking the large M limit of ( ^.3| ) 
Sl = liniM^oo jY^og Zl^m- Using the explicit form ( |4.18D , the largest eigenvalues are 
obtained by having the maximum number of terms {1 + zj) that are larger than one in 
module. 

We have represented in Fig. 21 the unit circle and its shifted image under z 1 + z. 
Only a portion of this image lies outside of the unit circle. Let us denote by z± = e^^*'^/^ 
the two intersections between the two circles of Fig.21. We will have a maximum number 
of terms (1 + Zj) with modulus larger than one for the value of ko such that 

^ < 1 < ^ (4.19) 
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Fig.21: The roots Zj of the Bethe Ansatz equation (|4.16D he in the complex 
plane on the unit circle (dashed circle). We have also represented the location 
of their shifted values (1 + Zj) (solid circle). The maximal eigenvalue among 
the values ( [4.18| ) is obtained by retaining only those Zj's whose shifted value 
lies outside of the unit circle (black dots) and having as many of these as 
possible. These correspond to Zj = e*^^ with 9j G [— 27r/3, 27r/3]. 

For simplicity, let us assume that L = 3n, in which case (|4.19|) gives ko = 2n. The largest 
eigenvalue reads then 

2n-l 

A2r^(0,l,...,2n- 1) = JJ (l + e2*^(^-^)/(3'^)) (4.20) 
The corresponding tiling entropy reads 

2n-l 

Sl = S3n=Yl Log(l + e2-(^-^)/(3-)) (4.21) 

4-5. Thermodynamic entropy and central charge 

To get the thermodynamic tiling entropy per tile ( |4.4| ), we must also take the large L 
limit of Sl- s = lim^^oo Sl/ L. We immediately get 

dxLog(l + 6^'''^) = > ^ ' 

-1/3 p>i 

The central charge may be obtained from the large L corrections to the thermodynamic 
entropy. From general principles, the large L limit is indeed expected to be described by a 
conformal field theory if the model is critical and translationally and rotationally invariant 
as it is presently. 

Moreover, the central charge c is extracted by computing the anomaly of free energy 
of the model on a cylinder of width L, that behaves like 
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C"' dx Log(l + e^--) = y sin ^ (4.22) 

J-i/3 ^ 3 



for large L, with a bulk free energy per tile / = — s as above and a central charge c. The 
non-conventional factor ^ is simply due to our deformation of the original rhombus of 
size L X M into a rectangle with size L x M-\/3/2 in the same units. 

The large L = Sn expansion of is easily derived using the following expansion, 
valid for any sufficiently differentiable function / 

2n — l 2n — l 2j + l 

-'-/'^ 2n-l ._ 2n-l STT 

J— 2 J~ 2 



J~ 2 



2T1-1 

2 



3(6n)"^ 3n n 

J~ 2 



E /"(5i;)+0(;;3) 
1 /-^/^ 1 

-— ^ / r(x)dx+o(^) 

6(6n)^ i_i/3 

f (1/3) - f (-1/3) 



with 



(4.24) 



6(6n)2 

Applying this to the function 

f{x) = Log(l + e^^^^^) (4.25) 
we finally get the expansion of 5"^: 



^ :(/'(l/3)-/'(-l/3)) = l (4.27) 



27rV3 

The above computation therefore confirms the result announced in Section 4.1 that the 
FPL{1) model is described in the continuum by a CFT with central charge c = 1. 



5. Effective field theory description via Coulomb gas 

5.1. Two-component Coulomb gas for FPL{2) 

Let us now see how to construct an effective continuum field theory for the FPL{n) 
model at large distances. We introduce a continuous description of the heights in the form 
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Fig.22: The triangular lattice T of mesh size a, dual to the hexagonal lattice 
on which loops are drawn, is the locus of the heights X allowed by the rules 
of Fig.20-(a). The ideal lattice X, corresponding to the values of /i = {X) 
where the free energy is minimal, is an hexagonal lattice of mesh size a\/?>. 
The triangular sub-lattice TZofX corresponds to ideal states associated with 
the same coloring. The mesh size of TZ is v^a. Passing from triangle M to 
triangle N by successive reflections along the edges brings the triangle back in 
the same orientation. The reciprocal lattice TZ* of the lattice 7?. is a triangular 
lattice of mesh size 2 / (3a) (here represented for a = 1/ \/3) . The smallest non- 
zero vector of the reciprocal lattice X* of X is the second smallest vector of 
TZ* , for instance that joining M to P in the figure. Its length is 2/{^/3a). 

of a two-dimensional two-component field h{r) = (/ii(r), /i2(r)) where r = (ri,r2) is the 
continuous version of the two coordinates of a node on the hexagonal lattice while h{r) is 
a locally averaged value of the two-dimensional height X around the point r. A precise 
definition of h is not necessary as only the symmetries of the problem will indeed matter. 



The approach presented here follows Ref.||3^ and is based on two-dimensional Coulomb 
gas techniques . 

A first step consists in identifying the so-called "ideal" configurations, i.e. those with a 
maximal entropy of excitations. We shall indeed assume that the large scale statistics of the 
problem is dominated by the fiuctuations around these ideal states. In the dual language of 
the planar folding of the triangular lattice, these ideal configurations correspond to maximal 
foldings, i.e. configurations where the lattice is completely folded onto a single triangle. 
These compactly folded states are indeed those with the largest number of elementary 
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excitations obtained by unfolding the six links bordering an elementary hexagon. To 
each ideal configuration, we may associate the corresponding averaged value of h, equal 
to the position of the center of mass of the triangle on which the lattice is folded. The 
corresponding allowed values of h form themselves a hexagonal lattice X of mesh size a / \/3 
if a is the mesh size of the folded triangular lattice (see Fig. 22). We finally write the 
effective free energy as 



f = J d\[ng {iVh^f + (V/is)') + Vih)] (5.1) 



with a Gaussian part describing the height fiuctuations and an effective "locking" potential 
V{h) having the periodicity of the ideal lattice X and enhancing the weight of the ideal 
states. The "roughness" g will be determined below in a self-consistent way. 

In terms of loops, each ideal state corresponds to a regular covering of the hexagonal 
lattice by means of elementary hexagonal loops of length six, all oriented in the same 
direction and centered on one of the three sub-lattices of the dual (tripartite) triangular 
lattice. There are six such tilings corresponding to the six permutations of the three colors 
ABC. The vertices of the ideal lattice X that correspond to the same covering, i.e. the 
same coloring, form a sub-lattice TZ of the ideal lattice X. This so-called "repeat" lattice 
7^ is a triangular lattice of mesh size ^/3a (see Fig. 22). Any local observable of the loop 
gas must therefore have the periodicity of the lattice TZ and may thus be written as 

$(r) = J2 $Ge^'"^-''("^ (5.2) 
Gen* 

where G runs over the vectors of the reciprocal lattice TZ* of TZ (defined by G.b = integer 
for any b in TZ). The lattice TZ* is a triangular lattice of mesh size 2/(3a) (see Fig. 22). 
On the other hand, the locking potential V{h) = 'YliOeTi* '^G^'^^^^'^^^^ has non-vanishing 
Fourier components vq only for G G X* , the sub-lattice of TZ* equal to the reciprocal lattice 
of the ideal lattice X. The value of g may now be obtained by assuming that V{h) is a 
marginal perturbation of the Gaussian free field, i.e. with dimension x = 2. 

Very generally, the dimension of the (electric) "vertex operator" Voir) = g2i7i-G./i(r)^ 
which describes the large distance behavior of the correlation function (VG(r)y_G(0)) ~ 
j,-2x(G) equal to |3| 

x{G) = (5.3) 
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The dimension x of V{h) corresponds to that of the most relevant electric operator (with 
the smallest x{G)) occurring in the Fourier decomposition of V{h), i.e. to the vector G 
of X* with minimal norm (vector MP in Fig. 22). The norm of this vector being equal to 
2/(-\/3a), we finally obtain 

*=3^ ^'■*'> 

Once g is fixed, we can obtain exponents describing the large distance behavior of 
correlations between coloring defects. These take the form of "magnetic operators" defined 
as follows. A coloring defect around a node of the hexagonal lattice corresponds to a 
dislocation-type defect in the height X after a complete turn around this node. This height 
defect reads AX = M = mi(A — S) +m2{A — C) for two relative integers mi and m2 (the 
vector M is thus a vector of the lattice TZ) . In order to have a configuration without defect 
at infinity, we must ensure magnetic neutrality, for instance by introducing another defect 
with magnetic charge —M at some other node. In this case, the variable X has a jump 
discontinuity equal to M along a line joining these two defects. The defect/anti-defect 
correlation then behaves at large separation r as r~'^^^^^^ where the magnetic dimension 
reads 

, , gM^ 1 [mV , , 

A first example of defect consists in having two unvisited sites at distance r from one 
another. This corresponds to a coloring defect characterized by M = 3A with norm 3a 
and with therefore x = 3/2. The fact that x < 2 shows that this is indeed a relevant 
perturbation and that the fully packed loop gas corresponds to an unstable phase with 
respect to such defects. The introduction of coloring defects is also useful to describe open 
lines in the loop gas. For instance, we may describe an open BC line linking two points 
at (an odd) distance r by introducing coloring defects at these two points, each with its 
three incident edges colored A,A,B. The corresponding value of M is M = 2 A + B, with 
norm -\/3a, which gives an exponent xi = 1/2. Very generally, one finds an exponent xe 
associated with the large distance behavior of correlations of the form r~^^^ for a set of i 
BC lines starting from the vicinity of a given point and ending at the vicinity of another 



at distance r, with the value 34 



X2k = -k 

I (5.6) 



according to the parity of i. 
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5.2. General FPL{n) and 0{n) cases 

Beyond the case n = 2, the FPL{n) model for —2 < n < 2 may be described 
by configurations of oriented self-avoiding fully packed loops, with weights e^"^*^/^, resp. 
g-i7re/6 pgj, right, resp. left turn of the loop. Summing over the two possible orientations 
of each loop results in a weight n — 2cos(7re) per loop. This holds only for contractible 
loops while loops escaping at infinity receive a weight 2 instead. To understand how to 
cure this problem, it is more convenient to define the model on an infinite cylinder, in 
which case the non-contractible loops wind once around the cylinder and receive a wrong 
weight 2 instead of n. This is corrected by introducing two electric operators with (two- 
dimensional) charges E and —E at both ends of the cylinder, which implies an additional 
energy term of the form 2i7rE.{h{+oo) — h{—oo)). Each non-contractible loop induces a 
discontinuity h{+oo) — h{—oo) of ±B (or ^C) and it is thus sufficient to choose E such 
that E • B = —E ■ C = e/2, i.e. E = e{B — C)/(3a^). Summing over the two orientations 
of the non-contractible loops now yield the correct weight n while contractible loops are 
unafi^ected. The introduction of the "background" charge E modifies the dimensions of 
the electric and magnetic operators into 

x{G,M) = ^G-{G-2E) + ^M^. (5.7) 

The value of g is again determined by the requirement that V be marginal. Choosing the 
determination < e < 1 of e, one finally gets 

One then finds an exponent Xi{n) associated to the correlation of £ BC lines given by 



1 - e 

1 - e ^2 
X2k-i{n) = — r— (fc^ - A; + 1) - 



(5.9) 



2 ' ' 2(l-e) 

Finally, the central charge of the model now reads c = 2 + 12x{E, 0), i.e.: 

Cfuiiy packed (?^) = 2 - ^Y^T^' ^ 2cos(7re) (5.10) 

To conclude this study, let us now see how the above formulas are modified in the 
case of the dense phase of the 0{n) model, which as we already saw for n = 2, consists 



44 



in allowing for the presence of unvisited nodes, in which case the heights become one- 
dimensional. The above analysis transcribed to this much simpler case leads finally to a 
central charge 

Cdensein) = 1 - Q- , n = 2cos(7re) (5.11) 

1 — e 

and to exponents 

Mn) = - 2(rri) (1 - ^^^o) (5.12) 

for any £. 




fully packed loops dense loops 




n 



Fig.23: Renormalization flow in the (n,w~^) plane. We have indicated an 
example of flow at constant n. The fully packed loop gas corresponds to an 
unstable flxed point at = 0. The dense loop gas corresponds to a stable 
flxed point. The central charge decreases by 1 between these two flxed points 

Cdense('^) — CfuHy packed ('^) 1- 

The lesson to draw from the above is that, at flxed weight n per loop, the fully packed 
loop gas has a central charge 1 more than the dense loop gas. This result is illustrated in 
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Fig. 23 where we display a diagram of renormalization group flows in the (n, plane, 
where u is the weight per visited site and —2 < n < 2 [^. All the flows follow lines of 
constant n and the problem is symmetric under u —u (since the number of sites visited 
by loops is even). The = line, which corresponds to the fully packed loop gas, is a line 
of unstable flxed points (see Fig. 23). The presence of defects (which, as we already saw, 
are relevant) drives the model away from this line towards the line of stable flxed points 
describing the dense phase of the 0{n) model. One also flnds another line of unstable flxed 
points at flnite values of = a/2 + a/2 — n corresponding to the transition points of the 
0(n) model [|0|. 
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PART C: LOOP MODELS ON RANDOM LATTICES 



We would like to generalize the models of 2D folding described in part A and those 
of fully packed loops described in part B to the case of random 2D lattices. By random 
2D lattices, we mean tessellations made of 2D polygonal rigid tiles glued together so as 
to form discrete surfaces with possible curvature defects concentrated at the nodes. We 
are interested in the phantom folding configurations of statistical ensembles of 2D lattices. 
These form discrete models for "fluid membranes", i.e. membranes without internal elastic 
skeleton and whose internal metric fluctuates 0. Beyond folding, we will also study fully 
packed loop models on random lattices. All these models are examples of a larger class of 
statistical models describing the coupling of "matter" degrees of freedom to the fluctuations 
of "space", used in the quantization of general relativity. More precisely, the present 2D 
models correspond to discrete realizations of the so-called two-dimensional quantum gravity 
(2DQG) 

In this part, we concentrate on the case of random triangulations or, dually, random 
trivalent graphs. Various results are presented, relying either on matrix integral techniques 
T9| , |39| - p] or more general effective 2DQG descriptions [Q. Exact predictions are tested 



against numerical enumerations . The lesson of this study is a clear distinction between 



ordinary 2DQG and so-called Eulerian 2DQG corresponding to a restriction of the set of 
dynamical lattices to Eulerian tessellations. 

6. Folding of random lattices 
6.1. Foldahility of triangulations 




Fig. 24: Example of triangulation with spherical topology. 
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In this part C, we shah consider only the simplest case of triangulations, for which 
all the tiles are equilateral triangles with edges of unit length. We will moreover restrict 
ourselves to planar triangulations, i.e. tessellations of a surface with the topology of the 
sphere (genus 0, no handles). An example of such a triangulation is given in Fig. 24. 

Given such a triangulation, we define as before a folding as a map from the trian- 
gulation into IR'^ such that its restriction to any elementary triangle is an isometry. In 
particular, each triangle is again mapped onto an equilateral triangle in M!^. 




f 



Fig.25: Two examples of planar triangulations: the planar representation 
(left) keeps track of the connectivity of the triangles only at the expense of 
deforming them. Note that the external face is also a triangle. The triangula- 
tion (a) is foldable into IR^ on a tetrahedron, but it is not foldable in the plane 
because its nodes are not tricolorable. The triangulation (b) is foldable into 
H on an octahedron and it is foldable in the plane as it is node-tricolorable 
(here with three colors e,f,g). An example of folding is indicated on the right 
in the planar representation by thickening the folded links. 

In the simplest case of planar folding (d = 2), a preliminary question concerns the 
very existence of such foldings. Indeed, not all triangulations may be folded into the plane. 
For instance, the triangulation on the top of Fig.25 (which may be embed in IR"^ so as to 
form a tetrahedron) cannot be folded into the plane. On the other hand, the triangulation 
at the bottom of Fig.25 (which may be embedded in JR^ so as to form an octahedron) can 
be folded into the plane. 

To answer this question, we note that, whenever a folded state exists, the positions in 
the plane of the nodes of the triangulation necessarily belong to a regular triangular lattice 
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of mesh size 1 in the plane H^. The triangular lattice is tripartite, i.e. its nodes may be 
colored with three colors (e,f,g) so that no two adjacent nodes be of the same color. This 
tricoloring is unique (up to a global permutation of the three colors). If a folded state 
exists, it therefore induces a tricoloring of the nodes of the triangulation we started from. 
This yields a necessary condition for the triangulation to be foldable in the plane, namely 
that it be node-tricolorable. Conversely, starting from a node-tricolorable triangulation, 
its tricoloring is unique (up to a global permutation of the colors). Such a triangulation 
may be folded on a single equilateral triangle by sending each node of a given color onto 
one of the three vertices of this triangle. Such folding corresponds to the complete folding 
of the triangulation from which other folded configurations may be obtained by partial 
unfolding. 

We obtain finally the following equivalent characterizations: 

1. The triangulation is foldable into the plane; 

2. The triangulation is node-tricolorable. 

In the case of a triangulation with spherical topology, we also have the following equivalent 
alternative characterizations ||43|| : 

3. Its faces are bicolorable (with distinct colors on adjacent triangles); 

4. Its edges may be oriented so that the boundary of each triangle receives a well-defined 
(clockwise or counterclockwise) orientation; 

5. The number of triangles around each node is even; 

6. The number of edges adjacent to each node is even. 

This last property justifies the denomination "Eulerian" for such triangulations as it 
ensures the existence on the triangulation of a closed Eulerian path (i.e. a path visiting all 
edges exactly once) . In other words, it ensures the possibility of drawing the triangulation 
by a single (closed) path without lifting the pen. 

To summarize, the planar foldable triangulations are the Eulerian triangulations, char- 
acterized by any of the above properties 1-6. 

6.2. Enumeration of foldable triangulations 



The enumeration of Eulerian triangulations was first carried out by W. Tutte 
in its dual version under the denomination of "bicubic maps". This enumeration may 
be generalized so as to keep track of the three colors, e.g. with different weights p, q 
and z for the nodes of color e, f or g respectively |^]. In other words, we enumerate 
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(a) (b) 

Fig.26: Equivalence between (a) a bicolored graph with arbitrary valences 
and (b) a tricolored triangulation. Prom (b) to (a), we simply erase all the 
nodes of the triangulation of a given color, say, g as well as all the links 
connected to them. Conversely, from (a) to (b), we note that the nodes 
around each face have alternating colors e and f. We then re-introduce at the 
center of each face a node of color g and connect it to all nodes around the 
face, thus creating triangles. 

random triangulations which are completely folded on a single triangle, keeping track of 
the numbers of nodes sent onto each of the three vertices of the triangle. 

This may be performed by use of matrix integral techniques. More precisely, we may 
compute the generating function Z{p, g, z; t; N) for possibly disconnected vertex-tricolored 
triangulations of arbitrary genus, with weights p, q, z for the three colors and a weight t 
per edge of color e-f. The parameter N governs the genus via a factor A^^ where x is the 
Euler characteristic of the triangulation at hand. The partition function is defined as 

Zip,q,z;t;N)= Yl p^-^^^g-^^^^.^-^^) ^ ^^1. (6.1) 

node — tricolored ' \ / I 

triangulations T 

where nej\g{T) denote the total numbers of vertices of color e,f,g, where A{T) denotes 
the total number of faces (twice the number of e-f edges), and where x(T) denotes the 
Euler characteristic of T. The factor |Aut(T)| is the order of the symmetry group of the 
tricolored triangulation T. 

The construction of a matrix model to represent Z{p, g, z; t; N) is based on the follow- 
ing simple remark: in a given tricolored triangulation T if we remove say all the vertices 
of color g and all edges connected to them, we end up with a bicolored graph, with uncon- 
strained vertex valencies (see Fig.26). The problem is therefore reduced to the enumeration 
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of bicolored graphs with a weight p (resp. q) per node of color e (resp. f), a weight z per 
face and a weight t per hnk. Such bicolored graphs are easily built out of the Feynman 
graphs of a two Hermitian matrix model, say Mg and M/, the index standing for the color. 
The colored m-valent vertices of the Feynman diagrams read 



Tr{Mr) = ( *:||:::: 

(6.2) 

for all m, and are to be connected via propagator edges with weight t/N while ensuring 
the alternation of colors e-f 

{{Ma)ij{Mb)ki) = {l-6ab)SjkSii^ = :=^::::: 

(6.3) 

where a,b = e or f . Let us introduce the corresponding matrix integral, but keep N fixed 
while the matrices are taken of size n x n, n possibly different from N. This gives the 
partition function 

Zn{p,q-t;N) = I [ dM,dMfe-^^'''^^^'^r,p,i,t) 

'pn{tj J-^ ) J 

V{Me,Mf;p,q,t) = pLog{l - M^) + qLog{l - Mf) + jM^Mf 

where the normalization factor (pn{t, N) ensures that Zn{0, 0;t; N) = 1. 
The Feynman graph expansion of the free energy reads 

Fn{p,q;t;N) = Log Z„(p, g; t; AT) 

V ^ (r) n f {V)+E{T) ^T{V{r)-E{T)) F{V) (6.5) 

^ |Aut(r)|^ ^ 



(6.4) 



bicolored conne-cted 
graphs r 



where we have denoted by na(r) the number of vertices of color a, V(r) = ne(r) + 
n/(r), -E'(r) the number of edges, F{T) the number of faces, and |Aut(r)| the order of the 
symmetry group of the bicolored graph F. Adding a central vertex of color g in the middle 
of each face of F, and connecting it to all the vertices around the face with edges will result 
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in a vertex-tricolored triangulation T. The number of such added vertices is nothing but 
rigiT) = F{r). Introducing 

. = ^ (6-6) 

we may rewrite 

Fr,{p,q;t;N) = F{p,q,z;t;N) (6.7) 

by use of the Euler relation 2 - 2h{T) = 2 - 2h{T) = V{T) - E{T) + F{T) and the fact 
that A(T) = 2E{r), as each edge of F gives rise to two triangles of T, one in each of 
the two faces adjacent to the edge. It is also a simple exercise to show that |Aut(T)| = 
|Aut(r)|. Hence computing Fn{p, q;t; N) through the integral formulation (|6.4|) will yield 
the generating function for compactly foldable triangulations. In particular, the generating 
function /(p, z; t) for planar tricolored triangulations with a weight q, z for the nodes 
of the three colors and t^^^ per edge is obtained by taking the N oo limit so as to select 
planar graphs only, and with n/N = z fixed, i.e. 

/(p, z; t) = lim -l^LogZ(p, q, t; zN, N) (6.8) 

An explicit calculation gives @ 

tdtf{p,q,z;t) = ^Wl^{l-U,-Uf-Ug) (6.9) 

where the functions Uf,, Uf and Ug are determined as formal power series of t by the 
equations 

U,{l-Uf-Ug) = pt 

Ufil-Ug-U,) = qt (6.10) 

Ug{l-Ue~Uf) = zt 

with the condition Ua = 0{t), a = e,f,g. This solution was alternatively recovered 
in Ref.^^ by solving a rectangular matrix model. Remarkably, the functions Ua/t are 
nothing but generating functions for rooted planar trees with tricolored vertices and a root 
vertex colored a. For instance, the first line of (|6.10| ) reads Ue/t = p/{l — t{Uf/t + Ug/t)) 
which is easily seen to generate all trees with root colored e weighted by p and arbitrary 
many descending subtree of color f or g attached to the root by inner edges weighted by 
t. The reason for this apparently mysterious coincidence between generating functions 
for node-tricolored triangulation and node-tricolored trees may be explained in a purely 



combinatorial way. Indeed, a series of works 46,47| has established bijections between 
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various classes of planar graphs and possibly decorated trees. In the present case, such a 
bijection exists which consist in cutting the triangulations so as to form tricolored trees 

mm- 

By eliminating Uf and Ug, we get the following 5th order equation for Ue 

Ul{l-Uef{l-2Ue + 2{p-q-z)t) = t'{{l-U,fp'-Ul{z-qf) (6.11) 

of which we should retain the unique solution with behavior ~ pt at small t. The values 
oiUf and Ug are obtained from Uf, by permuting the weights p, g, z. 




Fig.27: The value of q, z) in the domain z < q < p. 

The case p = q = z is much simpler as we may write Ue = Uf = Ug = U , where U is 
the solution of the quadratic equation U{1 — 2U) = zt, namely 

1 



U = -(1 - Vl-Szt) (6.12) 
The generating function f{z,z,z;t) therefore reads 

tdj{z,z,z;t) = ^{1-3U) = ^ (^(^2-^^ (l-v/r^8^) + l-6^t^ (6.13) 

from which we deduce the number Ng of Eulerian triangulations with 3s edges (i.e. 2s 
faces) and a marked oriented edge. 

2 s!(s + 2)! 2v^sV2 ^ ' 
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in agreement with Tutte's result [^. This result is to be compared to the case of arbitrary 
rooted planar triangulations with 2s faces, which behaves as (12v^)'^/s^/^ [|39|J4^ . 



We deduce from the solution ( |6.11|) that with arbitrary positive weights p, q, z, the 
free energy /s(p, z) of tricolored triangulations with 3s edges behaves generically for large 

s as 

Mp, q, z) ~ C(p, ^) MiilflL (6.15) 

S2 

where z) is the value of t at which the function becomes singular. By homogene- 

ity, we have ti,{p,q,z) = {l/p)T{z/p,q/p) and, by symmetry, we can restrict ourselves to 
the domain where, for instance, z < q < p. Figure 27 represents q, z) in this domain 
and gives some explicit formulas on the boundary of the domain. The value 7/2 of the 
exponent in formula (|6.15|) is a general feature of so-called "pure gravity" models that de- 
scribe the universality class of fluctuating 2D-space without matter. It is a universal result, 
common to all planar graph enumeration problems at generic values of their parameters, 
i.e. without fine-tuning. This is a particular application of the so-called KPZ relations 



described in detail in Section 7 below, for the coupling to gravity of a trivial conformal 
matter theory with central charge c = 0. 

The above enumeration of Eulerian triangulations solves the problem of compact fold- 
ing, i.e. enumerating configurations of random triangulations completely folded onto a single 
triangle. The Eulerian nature of a triangulation allows to orient all its links in a coherent 
way, which in turn allows for a proper definition of the link variables in the form of tangent 
unit vectors with a vanishing sum around each face. In order to describe the most general 
foldings of the triangulation, we simply have to assign colors to the links representing the 
orientation of the corresponding tangent vector in the folded configuration. To obtain 
all folded states, we have to consider all the tricolorings of these links with three colors 
A, B or C, distinct around each triangle. In other words, the folding of random planar 
triangulations is equivalent to the simultaneous tricoloring of the nodes and of the links. 
As opposed to the tricoloring of the nodes which is essentially unique, the tricoloring of 
the links is a source of entropy. While the problem of tricoloring of the links only was 
solved p0| , |51[1 , the simultaneous tricoloring of both the links and the nodes is still an open 



question. 
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7. Statistical models coupled to 2D Quantum Gravity 

The problem of folding of random triangulations studied in Section 6 above is a par- 
ticular case of statistical model coupled to 2D Quantum Gravity (2DQG). This simply 
means that we replace the underlying regular lattice of some ordinary 2D statistical model 
by a somewhat arbitrary tessellation of the sphere or of some higher genus surface. The 
tessellation becomes therefore part of the configuration to be summed over. Such models 
were introduced as discrete descriptions of 2D quantum gravity, where "matter" models 
are coupled to the quantum fluctuations of the underlying 2D "space" , represented in the 
discrete by the tessellations (see e.g. Ref. for a review and more references). To each 



tessellation G, we associate a statistical weight directly borrowed from the Einstein action 
in 2D, namely a weight A^x(©)(^^(0) where x ^-nd A are respectively the Euler characteristic 
and area of the tessellation (respectively measured via x = #faces— ^^edges+^^vertices and 
A = ^faces or vertices), and where N and g are the discrete counterparts of the Newton 
constant and cosmological constant. For each tessellation G together with a matter con- 
figuration, realized for instance via a set of spins {a} on G, we also have a weight e~^^^'^^^ 
for some energy functional E. We also make the standard choice of dividing the resulting 
weight by the order |Aut(G, {cr})| of the automorphism group of the tesselation together 
with its spin configuration. This leads for instance to the discretized partition function of 
any 2D statistical lattice model coupled to 2D quantum gravity 

Z= Yl iV^^''^(7^(e) J]e-^^«'^>V|Aut(G,{a})| (7.1) 

tessel. {cr} 

Note that the parameter allows to isolate the contributions of tessellations of fixed genus. 
In the following we will be mainly interested in the genus zero contributions, obtained by 
letting oo. Note also that the "free energy" F = LogZ selects only the connected 

tessellations in the sum ( |7.1D . 

Like in the fixed lattice case, we are interested in the thermodynamic limit of the 
system, in which say the average area or some related quantity diverges, ensuring that the 
dominant contributions to Z come from large tessellations. This is guaranteed in general by 
the existence of a critical value Qc of the cosmological constant g at which such divergences 
take place. This value is a priori a function of the type of random lattices we sum over as 
well as of the various matter parameters. We may now attain interesting critical points by 
also letting the matter parameters approach critical values, a priori distinct from those on 
fixed lattices. The result is well described by the coupling of the corresponding CFT's with 
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quantum fluctuations of space, namely by letting the metric of the underlying 2D space 
fluctuate. Such fluctuations may be represented in the conformal gauge by yet another field 
theory, the Liouville field theory, which is coupled to the matter CFT. This field theoretical 
setting allows for a complete understanding of the various critical exponents occurring at 
these critical points [^2|. For instance, one defines the (genus zero) string susceptibility 
exponent --ystr as the exponent associated to the cosmological constant singularity, namely 
by writing the singularity of the free energy as 

i^|singOc((7c-^7)'-^^"^ (7.2) 

In the case of coupling of a matter theory with central charge c to 2DQG, one has the 
exact relation 



c-1- v/(l-c)(25-c) 

Istr = IstAc) = (7.3) 

In the case of "pure gravity", namely when the matter is trivial and has c = 0, we get 
Istr = ~l/2, while for the critical Ising model with c = 1/2 we have 'jstr = —1/3. 

Upon coupling to gravity the spinless operators of the CFT {^h{z,z)) get "dressed" 
by gravity ($/i = ^a) and acquire dressed dimensions A, given similarly by 



V25-C- ^ ' ' 

As opposed to the fixed lattice case, where conformal dimensions govern the fall-off of 
correlation functions of operators with distance, the dressed operators of quantum gravity 
do not involve distances, as their position is integrated over the surfaces, but rather only 
involve changes of area at fixed genus. More precisely the general genus zero correlators 
behave in the vicinity of [p2[ as 



(vE'A,^'A,...^'A,) ~ {9c - gf ^^*'-+^i<-^^^' (7.5) 

These results may be easily translated into the large (but fixed) area A behavior of 
the various thermodynamic quantities, upon performing a Laplace transform, which selects 
the coefficient of in the various expansions. Let Fa denote the partition function for 
connected tessellations of genus zero and area A, we have 
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while if {■■■)a denotes any genus zero correlator at fixed area, we have 

(v&a,^a....^a.)a ~ 3 ^ (7.7) 

In Ref ■ |[42|| , all these formulas were also generalized to higher genus as well. Note finally 
that Eqs.( [7.3|) and ( |7^ ) are valid only as long as c < 1. This corresponds to the famous 
"c = 1 barrier" beyond which the gravitational system degenerates into infinitely branched 
structures (branched polymer phase of 2DQG). 

As an illustration of Eq.( |7.6| ), recall that the number of Eulerian triangulations of the 
sphere with A triangular tiles and with a marked edge reads from Eq.( |6.14| ) 

N 3,^/, {A)\ [2 8^/^ 

with A = 2s the number of triangles. Noting that the rooting simply amounts to A^^ oc 
AFa, hence the asymptotics (|7.8[ ) correspond to gc = l/VS and 'jstr = —1/2, hence c = 
according to (|7.3| ). This is one of the various examples where we attain the universality 
class of pure gravity, namely by summing over some bare tessellations without matter on 
them. 

The asymptotic enumeration for most combinatorial problems involving planar (or 
more generally fixed genus) graphs is encoded in Eqs. (|7.6|) and ( [7.7|) or their higher genus 
generalizations. This allows to predict the corresponding configurational exponents pro- 
vided one is able to identify the central charge c of the underlying GET. This latter step 
however may prove to be quite involved. In fact, a lesson of the forthcoming Section is 
that the naive application of these formulas may lead to wrong results for statistical mod- 
els whose definition strongly relies on the structure of the underlying (whether fixed or 
random) lattice. This is precisely the case for the problems of folding or of fully packed 
loops that we are interested in. 



8. One-flavor fully packed loops 

8.1. Fully packed loops on random trivalent graphs 

We now come to the random version of the FPL{n) model, i.e. to a gas of fully 
packed loops on random planar trivalent graphs (dual to random triangulations). The most 
natural version consists in summing over all planar trivalent graphs. A typical, connected 
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Fig.28: An example of configuration of fully packed loops on a random planar 
trivalent graph. 



configuration of this model is represented in Fig.28. It is made of a set of disjoint or nested 
loops, each with a weight n, linked together by non-crossing arches, corresponding to the 
links not visited by loops. 

The partition function Zg of the model, which counts the number of configurations 
with 2s nodes behaves as 

Zs-C^ (8.1) 

where C, R and 7 depend on n. The value of 7 was obtained exactly by matrix integral 
techniques (see Ref. [|2l) with the result 



7 = n = 2cos(7re), 0<e<l (8.2) 

1 — e 

We may also consider the partition sum counting FPL{n) configurations in the 

presence of two marked points connected by / open lines. The latter is characterized by a 
configurational exponent ai{n) through ~ / s"'-^'^\ with the result [53] 



ai{n) = l+^- (8.3) 

independently of n. 

As an illustration, a particularly interesting case is that of n for which configu- 
rations have only one loop. This case therefore describes the problem of enumeration of 
Hamiltonian cycles, i.e. closed loops passing through all the nodes of the random trivalent 
graph [|10|- If we open the loop at any of the 2s visited links and stretch it into an infi- 



nite oriented line as in Fig. 29, we simply obtain a configuration with 2s nodes on the line 
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Fig.29: An example of oriented Hamiltonian cycle on a planar random triva- 
lent graph (left). By cutting the link and stretching the loop into an oriented 
infinite straight line, one obtains a system of arches on top and below the line 
(right). 



connected by pairs via a set of, say k non-intersecting arches in the upper half plane and 
s — k in the lower half plane. We deduce that 



2sZ, 



s 

E 

A;=0 



2s\ 

2ik j CgCg-i-i 



4 16" 



(8.4) 



where the Catalan number Cg = (2s)!/(s!(s + 1)!) counts the number of possible configura- 
tions of s non-intersecting arches. One thus finds an exponent 7 = — 1, which is compatible 
with Eq.(|3) forn = (e = 1/2). 




1=1 




Fig.30: Examples of configurations with two marked points connected by / 
open lines, for / = 1 (top) and / = 2 (bottom) 
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We may similarly count configurations with two marked points connected by / open 
lines. For instance, the case I = 1 corresponds to having a segment dressed by arches 
(see Fig. 30). The number of such configurations with 2s trivalent nodes is clearly given 
by 2^*Cs ~ in agreement with Eq.( p.3|) for / = 1. Similarly, for / = 2, one 

describes configurations of a Hamiltonian cycle with two marked links (see Fig. 30), in 
number As'^Zg ~ in agreement with Eq. (|8.3| ) for I = 2. 

By inverting the relation (|7.3|) , we see that the value ( |8.2|) of 7 is compatible with a 
central charge c = Cdense{n) = 1 — 6e^/(l — e). Similarly, by using Eg. ([f77|) with /c = 2, 
i.e. ai{n) = 1 — 7 + 2Ai, and by inverting the relation (|7.4| ) , one finds that the result ( |8.3| ) 
for ai{n) is compatible with a conformal dimension 2h = xi{n) as given by Eg. ( ^.12]) . To 
conclude, the model FPL{n) defined on arbitrary random trivalent graphs corresponds to 
the coupling to gravity of a conformal theory of central charge Cdensein) as if the loops 
were not fully packed, a result further confirmed by the identification of the spectrum 
of dimensions xi{'n). Naively, we would have expected a central charge egual to that, 
Cfuiiy packed (^) 5 of the FPL{n) model on a regular lattice. In other words, we do not 
observe here the phenomenon c — c + 1 which we found for the regular lattice when 
passing from dense loops to fully packed loops. 

The explanation for this clearly comes from the strong coupling of the matter model 
to the symmetries of its underlying space. For instance, fully packed loops on the regular 
lattice all have even length, a necessary condition for viewing them as the BC loops of 
some ABC tricoloring. We can therefore consider a slightly constrained model in which we 
impose that loops be of even length on the random lattice. It was shown however ||5^ that 



the 0{n) model with loops of even length on random graphs is in the universality class of 
the 0{n/2) model. This does not lead to the desired c — c + 1 phenomenon. As we shall 
see just below, we need to apply a much more drastic constraint by reducing the class of 
random graphs itself in order to recover the desired central charge CfuUy packed (^) of the 
regular case. 

8.2. Fully packed loops on random trivalent bipartite graphs 

On the regular lattice, the central charge c = 2 of the FPL{2) model as opposed 
to c = 1 for the 0(2) model had a simple geometrical interpretation, as describing the 2 
degrees of freedom of the 2d positions of nodes for the eguivalent 2d folding problem of 
the triangular lattice. On the other hand, we also observed that this 2d folding picture 
could be extended to the case of random triangulations provided we restricted their class 
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to the Eulerian ones. Indeed the latter are precisely the triangulations foldable onto the 
plane. In the dual language of trivalent graphs, this condition amounts to requiring that 
the graph be bipartite, i.e. may be node bicolored. Beside the folding interpretation, we 
also noticed in the Coulomb gas approach of Section 5.1 that the rules for defining a two- 
dimensional height use explicitly the bipartite nature of the hexagonal lattice (see Fig. 20). 
We therefore expect this bipartite nature to be also crucial in the random case. 

In this Section, we concentrate on the two cases n = and n = 1 for which we show 
that the FPL{n) model, when defined on trivalent node-bicolored graphs, does indeed 
have the "increased" central charge CfuHy packed (^) = Cdense(^) + 1 as given by Eq.( |5.1Cl| ). 
In the case n = 0, this translates into a remarkable irrational critical exponent for a very 
simple apparently harmless combinatorial problem. 

For n = 0, the FPL{0) model describes the configurations of a Hamiltonian cycle, i.e. 
a single loop visiting all the nodes of the random bicolored trivalent graph (or equivalently 
visiting all faces of the Eulerian triangulation) . We would like to compute the number 
of such configurations with 2s trivalent nodes. 




Fig.31: An example of configuration of arches around a line obtained by 
cutting the Hamiltonian cycle. The white and black nodes alternate along 
the line. Each arch connects a white node to a black one. 



Once again we may open the Hamiltonian cycle at any of its 2 s links so as to build an 
oriented line containing all the original nodes, together with a system of arches connected 
them by pairs either in the upper- or lower-half plane delimited by the line (see Fig.31). 
Comparing this situation with that of previous Section, we now have two additional con- 
straints: 

(i) the colors of the nodes along the line alternate between black and white; 

(ii) each arch connects nodes of opposite colors. 

This combinatorial problem, although extremely simply stated, is still open. We can 
however estimate the configuration exponent of the problem with reasonable accuracy 



by direct enumeration for small enough sizes s. A first approach, explained in Ref.||Tl 
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consists in first building the two systems of arclies on top and below the line and to count 



the number of ways to interlock them so that the black and white nodes alternate. A 
second approach consists in using a transfer matrix to generate the systems of arches from 
left to right. A similar technique will be discussed in detail in part D devoted to meanders. 
The first values of sZf for s = 1, 2, ■ ■ ■ , 22 are given in the table below: 



1 


2 


2 


8 


3 


40 


4 


228 


5 


1424 


6 


9520 


7 


67064 


8 


492292 


9 


3735112 


10 


29114128 


11 


232077344 


12 


1885195276 


13 


15562235264 


14 


130263211680 


15 


1103650297320 


16 


9450760284100 


17 


81696139565864 


18 


712188311673280 


19 


6255662512111248 


20 


55324571848957688 


21 


492328039660580784 


22 


4406003100524940624 



We may easily extract from these data an estimate for the exponent 7 (see Fig. 32), 
namely 7 = —0.77(1). This value differs clearly from the value 7 = — 1 obtained in previous 
Section for non-Eulerian triangulations. Moreover, this new value is now compatible with 
a central charge c = CfuHy packed (0) = —1, for which, according to ([7.3|), 7 reads 

^ ^ _ 1 + VU 0.76759- ■• (8.5) 

6 



More precisely, by inverting the relation (|7.3D, our estimate for 7 yields a central charge 



c = — 1 ± 0.05. This leaves not much doubt on the fact that c = Cfuiiy packed(O) = —1 as 
expected. As a consequence, we predict that our simple combinatorial problem displays a 
quite remarkable irrational configuration exponent, given by Eq.( |8.5|) . 

We now turn to the case of the FPL(l) model on node-bicolored trivalent random 



graphs, discussed in Ref.[|T2|. For n = 1, and only at this value, the connectivity of 
the loops plays no role and allows for local transformations. Let us orient all the links 
toward their adjacent black node. By contracting into a single tetravalent node each pair 
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Fig.32: Estimates of 7 extracted from the exact enumeration for sZ^ (see 
table) in the Eulerian (o) case and from 2sZs = CsCs+i in the non Eulerian 
case (x). The x-axis coordinate corresponds to the minimal value of s used 
in the determination. The different curves correspond to successive iterations 
of the convergence algorithm used (see Ref. [|lll) . The value of 7, equal to 



T in the non Eulerian case, is equal to ~ —0.77(1) in the Eulerian case. 

of trivalent nodes separated by an unvisited link, we end up with a tetravalent graph 
with oriented links. Moreover the orientations obey the so-called ice-rule of the six-vertex 
model, namely that each node have exactly two ingoing and two outgoing arrows (see 
Fig. 33). More precisely, the two ingoing arrows must be consecutive around the node, 
which corresponds to the particular point of the six-vertex model where one of the weights 
[ws) is zero (see Fig. 33). The six- vertex model was solved by use of matrix integrals in 
Ref.fS^, where it was shown on one hand that the ws = point corresponds to a critical 



point, and on the other hand that the latter is described by the coupling to gravity of 
a particular CFT with central charge c = 1. This corresponds precisely to the expected 

result C = Cfully packed (n = 1) = 1. 

To conclude this part, the above analysis of the cases n = and n = 1 suggests that the 
specific universality class of fully packed loops observed on the regular lattice is preserved 
provided we impose the bicolorability of the underlying trivalent graph, i.e. the Eulerian 
nature of the dual random triangulation. Summing over such triangulations will be called 
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Fig.33: Mapping of the FPL{1) model on node-bicolored trivalent graphs 
onto a particular point of the gravitational six-vertex model on random 
tetravalent graphs. Top: orienting all links toward black nodes and con- 
tracting unvisited links, we obtain tetravalent nodes with two consecutive 
ingoing and two consecutive outgoing arrows. Bottom: the six local environ- 
ments of the "six-vertex model" and their respective weights on the regular 
square lattice. Once defined on a random tetravalent graph, the first two 
pairs of environments are indistinguishable, hence we must set wi = W2- The 
constraint of having two consecutive ingoing arrows imposes ws = 0. 

Eulerian gravity for obvious reasons. Now we may define two possible "gravitational" 
models of fully packed loops on random trivalent graphs. 

• We may sum over arbitrary trivalent graphs. As the dual triangulation will in 
general not be Eulerian, it will not be foldable, and the extra degree of freedom will be 
lost. The string susceptibility of the corresponding gravitational model will be computed 
using Eq.( |7.3|) with the dense central charge 



ordinary gravity : c — 1 — 6 



1 



n 



2 cos 7re 



.6) 



• We may sum over bipartite trivalent graphs, whose dual triangulation is automat- 
ically foldable, thus preserving the height variable in IR^. The string susceptibility must 
be computed using Eq. (|7.3| ) with the fully packed central charge 

„2 



Eulerian gravity : c = 2 — 6 



1-e 



n = 2 cos 7re 



(8.7) 



In the case n = 2 describing the 2d folding of (foldable) random triangulations, we ex- 
pect a GET with central charge CfuUy packed (2) = 2 whose coupling to gravity would lead us 
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beyond the above-mentioned c = 1 barrier above which the relation ( [7 .31) no longer applies. 
On heuristics grounds, we expect that the statistics of folding should be dominated by con- 
figurations where the folded surface degenerates into a highly branched one- dimensional 
structure. A confirmation of this image would require solving the simultaneous tricoloring 
of both nodes and links of random triangulations, an open question. 
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PART D: MEANDERS 



This part is devoted to the study of meanders, also equivalent to self-avoiding fold- 
ing configurations of a one-dimensional lattice. Beyond exact results for meander-related 
problems, we present a random-lattice fully packed loop description of the problem leading 
to precise predictions for the asymptotics of various meandric numbers. These predictions 
are tested numerically with remarkable agreement. A lesson of this study is the equivalence 
between the ID self-avoiding folding and the 2D phantom folding of random foldable but 
non-unfoldable quadrangulations. 

9. ID self-avoiding folding: Meanders 

The study of previous Sections was limited to phantom folding problems, with lattices 
made of interpenetrable cells. In particular, wc concentrated on statistical models of folded 
states without any reference to an actual folding process in a possibly higher dimension. 
The question of self- avoidance, where we now prevent cells from interpenetrating one- 
another is extremely difficult. As we shall see now, even in the simplest ID case, the 
problem is highly non-trivial and turns out to belong to the same class as the so-called 
meander problem, a notoriously difficult subject. 

9.1. The meander problem 




Fig.34: The M3 = 8 configurations of meanders for 2n = 6 bridges. 

The problem of meanders is one of the combinatorial problems which, although very 
simply stated, still resist any attempt to an exact enumeration, would it be only for 
the asymptotics. The meander problem may be stated as follows: find the number M„ 
of topologically inequivalent configurations of a closed non-intersecting circuit crossing a 
river through 2n bridges. In this formulation, the river is assumed to be an infinite oriented 
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line. Note that, by closing the river into a loop and by opening and deforming the circuit 
into a line, the circuit and the river play symmetric roles. The denomination "meanders" 
refers rather to the dual picture where the river meanders around the circuit but, for 
historical reasons, we shall work within the first formulation. For illustration, the figure 
34 represents the M3 = 8 configurations of meanders with 6 bridges. In the following, we 
shall be mainly interested in the asymptotics of for large numbers of bridges, as well 
as various generalizations. 

9.2. Meanders as a ID self-avoiding folding problem 




circuit 



Fig.35: Correspondence between (left) the folding of a dosed strip of 2n 
stamps and (right) a meander with 2n bridges. For clarity, the stamps are 
represented with different lengths and the folding is slightly undone. 

The problem of meanders is related to that of the self-avoiding folding of a closed 
chain of segments of unit length. Imagine for instance a strip of 2n post-stamps attached 
so as to form a loop. We want to count the distinct ways of maximally folding this strip 
on top of a single distinguished stamp. The main difficulty of the problem comes from 
the avoidance of the stamps which cannot interpenetrate one-another. Such a folding is 
represented on the left in Fig.35. We may now imagine piercing the stamps with a needle 
and a piece of thread and then knotting the thread into a loop (see Fig.35). Let us then 
open the strip at the level of the marked stamp and unfold it into a straight line (river), 
while the thread (circuit) meanders around it (see Fig.35). This transformation is clearly 
a one-to-one mapping between the ID self-avoiding foldings of a closed chain of length 2n 
on top of one of its segments and meanders with 2n bridges. 
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semi-meander 



Fig.36: The four ways (a) of folding a strip of 3 stamps (here viewed from 
the side) on a single stamp. The transformation (b) of a compact folding of 
a strip of n — 1 stamps (onto a single stamp) into a semi-meander with n 
bridges. The semi-meanders correspond to the case of a semi-infinite river 
with a source (•) around which the circuit may wind. 

The same correspondence applies to the folding of an open chain of, say n — 1 segments 
attached to a fixed support (see Fig.36-(a)). In this case, the transformation leads to what 
is called a semi-meander with n bridges (in correspondence with the n — 1 segments and 
the support), i.e. a configuration of a closed non-intersecting circuit crossing a semi-infinite 
river (with a source) through n bridges. Note that in this case, the circuit may wind freely 
around the source. Both meanders and semi-meanders will be studied below in a unified 
framework. 



9. 3. A brief history of meanders 

The meander problem is quite old. It is mentioned as early as 1891 under the 



name of "probleme des timbres-poste" (post-stamp problem) by E. Lucas [^. Later on, 
Sainte-Lagiie devotes a chapter of his book "Avec des nombres et des lignes: recreations 
mathematiques" (with numbers and lines: mathematical entertainments, 1937) [^. Sev- 
eral combinatorial approaches of this problem are discussed in Refs. ||58| - |60[] . More recently, 
the problem re-appeared with the work of Arnol'd |^l| in connection with the 16th Hilbert 
problem (the enumeration of ovals of algebraic planar curves). The modern formulation of 
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the problem, as just presented above, is due to Lando and Zvonkin |6^, who also intro- 
duced the name "meander" . Let us finally mention that the meander problem has many 



facets in relation with: mathematical problems such as the classification of 3-surfaces |]63 
computer science problems such as the study of planar permutations p^, or even artistic 
questions such as the description of mazes in Roman mosaics ||65|| . 

Despite a number of attempts, the problem is still unsolved to this day. The most 
important result of this part D is the prediction of a number of configurational exponents, 
all irrational, for meanders and such. 



10. Solvable cases 



A number of exact results are known for variants of the meander problem. These 
results rely on various descriptions using either direct combinatorics via arch statistics fl^ , 



explicit evaluations via matrix models ||6^] or algebraic formulations via the Temper ley- 



Lieb algebra [14 . 



10.1. Generalization: multi-circuit meanders 
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In the following, we consider the generalized case of meanders with k circuits still made 
of a single river but with a number k of non-intersecting, possibly interlocking, circuits. 

(k) 

We shall denote by the number of topologically inequivalent configurations of k non- 
intersecting circuits which cross a river through a total number of 2n bridges, with clearly 
k <n (we require that each circuit crosses the river). More simply, we shall refer to these 
configurations as meanders with 2n bridges and k circuits (see Fig. 37). Similarly, we shall 

(k) 

denote by SMn the number of semi-meanders with n bridges and k circuits corresponding 
to the case where the river is a semi-infinite line with a source around which the k circuits 
may wind (see Fig. 37). Here again, one necessarily has k < n. With these definitions, the 
original problems of meanders and semi-meanders correspond to /c = 1. 

Instead of working at constant k, it is easier to let k vary and to introduce a weight 
q per circuit. One then defines the partition functions 



Mn(g) = E^'^n'^ 
k=0 

n 



k=0 

The original cases of meanders and semi-meanders can be recovered by considering the 
limit 9 ^ of Mn{q)/q (resp. SMn{q)/q). 

For large n, we expect asymptotic behaviors of the form: 

MM-C{q)^ 

I Z (10-2) 

where a{q) and a{q) are the configurational exponents. In the case of semi-meanders, we 
can also define the winding w as the depth of the source, namely the minimal number of 
circuit crossings in a path from the source to infinity (see Fig. 37). We may then consider 
the average value of w for configurations with n bridges and a weight q per circuit. We 
expect the following behavior at large n: 



(«;)„(9)~n^(«) (10.3) 
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Fig.38: By opening the river, we transform any semi-meander with n bridges 
into a meander with 2n bridges. The opened circuits are completed by con- 
necting diametrically opposite bridges via nested arches forming a rainbow, 
thus keeping track of their connectivity. The winding (here w = 3) corre- 
sponds to the number of arches above the middle point (•). 

with a winding exponent < v{q) < 1. It is clear that n and w have the same parity and 
that meanders correspond to semi-meanders with 2n bridges and with a winding = 0. 

Note conversely that a semi-meander with n bridges may be viewed as a particular 
meander with 2n bridges by "opening" the river (see Fig.38), hence splitting each bridge 
in two. The connectivity of the loop is preserved by connecting the newly formed pairs of 
bridges via a "rainbow" of nested arches. 

Two situations are a priori possible: 

(i) ^{q) < 1 and R{q) = R{q). This case corresponds to a situation where the winding 
becomes negligible for large n and the entropy per bridge for meanders and semi- 
meanders is thus the same; 

or: 

(ii) h'{q) = 1 and R{q) < R{q)- This case corresponds to a situation where the possibility 
of winding generates an extra configurational thermodynamic entropy per bridge. 

We shall see in next Section that both situations (i) and (ii) actually occur, depending on 
the value of q. 

10.2. Combinatorial solutions at q — oo, 1 and —1: exact enumeration via arch statistics 

The list of exactly known enumerations of multi-circuit meanders reduces to the three 
cases q —>■ oo, q = 1 and g = — 1. In all cases, the problem is reduced to that of arch 
statistics. A multi-circuit meander with 2n bridges may indeed be seen as the juxtaposition 
of two arbitrary systems of n arches, one on each side of the river, and connected at the 
bridges. Similarly, by opening the river (see Fig.38), any multi-circuit semi-meander is 
the juxtaposition of a system of n arches on one side and the particular system made of a 
rainbow of n nested arches on the other side. 
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The number of possible configurations of n non-intersecting arches is given by the 
celebrated Catalan numbers: 

= (10.4) 

n\ (n -(- Ij! 

In the case Q = 1 where we ignore the number of circuits, we deduce that: 

1 42"^ 



M^{q = 1) = (c„)2 
SM^{q = 1) = c„ ~ 



r\j 

3 



I 4n (10.5) 



7r 77,1 



from which we get -R(l) = -R(l) = 4, = 3 and = 3/2. The winding number w 
corresponds to the number of arches passing above the source of the river (middle-point 
in Fig. 38). Its average value is easily evaluated, with the asymptotic behavior: 

(w),(q = 1) ~ (10.6) 



i.e. fi^q = 1) = 1/2. This is a particular instance of the situation (i) of previous Section, 
in which the winding is negligible. 

In the case q — * 00, we have to maximize the number k of circuits for fixed n, i.e. 
take k = n. In the case of meanders, this amounts to requiring that the systems of arches 
above and below the river be identical up to reflection. Similarly, in the case of semi- 
meanders, the system of arches above the (opened) river must be made of n nested arches. 
By re-closing the river, this selects a unique configuration made of n nested circles winding 
around the source. This gives: 

1 (27^)2- 



MM ~ Criq 



n 



/yf n-^/^ (10.7) 
SMM q^ 

and thus R{q) — * 2^^, R{q) q, a{q) 3/2, a{q) 0. One has clearly {w)n{q) n 
as the unique semi-meander with n circuits has winding n, and thus u{q) = 1. This now 
corresponds to the situation (ii) described in previous Section. We shall see later how this 
q ^ 00 result may be taken as the starting point of a systematic 1/q expansion. 

Finally, a last solvable situation concerns the case q = —1 for which one shows that 

, . / n even 

n even 



SM^{q = -l] 



—Cp n = 2p + 1 



This result may be proved by use of a simple involution for arches [|T^, or by a more 
technical supersymmetric matrix model approach |^ . 
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10.3. Multi-river, multi-circuit meanders: asymptotic enumeration via matrix models 

We now address a generalization of the meander problem in which we allow for arbi- 
trarily many rivers forming a set of nested or disjoint loops, crossed by arbitrarily many 
non-intersecting circuits, each with a weight q. We shall denote by M^(l, q) the generating 
function of the corresponding connected configurations with one of the rivers opened into 
a line. 

Following the work of Ref. |6^, we may use a "black and white" matrix model to 
evaluate the exact asymptotics of Mnil^q). We proceed by computing the generating 
function Zi^q{N;x) for possibly disconnected configurations of arbitrary topology, with a 
weight X per bridge and the standard gravitational weight A^^, x the Euler characteristics 
of the corresponding graph. For integer q, this function is given by the multi-matrix 
integral 

Zi^q{N;x)= dWdBi---dBqe V a=i a=i J (10.9) 

where the integral is over Q -f- 1, N x N Hermitian matrices and normalized so that 
Zi^q{N;0) = 1. The quantity M^il^q) is recovered in the planar limit ^ 00 as the 
coefficient of x^"^ of the quantity dxljOg{Zi^q{N; x))/N'^ . Indeed, the diagrammatic expan- 
sion of ( |10.9| ) involves connecting black or white half-edges representing matrix elements 
of Ba or W into a closed graph via the propagators 

white edges : {WijWu) = j^^u^jk = 
black edges : {{Ba)ij{Bh)ki) = j^Sa,bSiiSjk 

and with simple intersection vertices 



(10.10) 



Ay 



Tr{WBaWB, 



— f • — ^ 

Ay 



(10.11) 

This results in graphs made of black loops intersecting white ones, and the sum over the 
index a produces a weight q per black loop. This allows to identify white loop with rivers 
and black ones with circuits. 
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Performing the Gaussian integration over all Ba matrices, we are then left with 

This is just a Gaussian average over one Hermitian matrix W. 

The critical singularity of the genus zero free energy / = limTv^oo Log(Z'i ^(A^; x))/N'^ ^ 
{x{q) — x)^""*" is found to lie at a critical value x{q) given by 

^(?) = ^ ■ t e, (10-13) 

2sm (7r|) 

where we have set q = 2cos(7re) (0 < e < 1), while the corresponding critical exponent 7 
reads 

7 = -T-^ (10.14) 

1 — e 



This exponent is precisely that expected from the KPZ formula (|7.3|) for the coupling to 
2DQG of the dense 0{q) model with central charge Cdense(q') = 1 — 6e^/ (1 — e) as given by 
Eq. ( ^.111) . This shows that meander problems may be viewed as the coupling to gravity of 
particular critical loop models, and that the universality class of the present case coincides 
with that of the dense 0{q) loop gas. 

The above results translate into the multi-river meander asymptotics 



with 

1 sin^(7T„ 



1 _o^m^(^f) 
x{q) 



(10.16) 

a{l,q) = 

1 — e 

In particular, we recover from these values the case of meanders with arbitrarily many 
rivers and one single circuit by taking = 0, e = 1/2, and i?(l, 0) = 4, a(l, 0) = 3. These 
values coincide with those of previous Section for one river and arbitrarily many circuits, 
as it should by river-circuit duality. We list a few of the values R{l,q) and a{l,q) for 
various fractions e in the table below. 
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Fig.39: The QniQ) matrix (here for n — 3, — 5) obtained by juxtaposing 
all pairs (^4, B) of systems of n arches. The matrix element Gn{Q)AB = g^*^"^'^^ 
encodes the number of circuits c{A, B) of the obtained meander. 



Beyond the particular enumerations above, another exact result concerns the compu- 
tation for any q of the meander determinant defined as follows. We have seen that any 
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meander with 2n bridges (and with an arbitrary number of circuits) may be viewed as the 
jiixtaposition of two systems of n arches, say A and B. By taking aU pairs of systems of n 
arches, we define a symmetric square matrix Qnil) of size Cn x Cn whose element ^n(<?)yis, 
indexed by the pair {A, B) , is equal to 



Qn{q) 



AB 



(10.17) 



where c{A, B) denotes the number of circuits of the meander obtained by jioxtaposing the 
arch systems A and B (see Fig. 39). We shall call "meander determinant" the determinant 
of the matrix Qn{(l)- It is remarkable that this determinant may be computed exactly, with 
the result 



i=l 



2n 
n — i 



2n 

21 . J + 

,n — z — 1 



2n 
n — i — 2 



(10.18) 



where the C/^'s are Chebyshev polynomials defined recursively as Uj+i{q) = qUj{q) — 
Uj-i{q), Uo{q) = 1, Ui{q) = q. Equivalently, we have Uj{2 cos{9)) = sm{{j + 1)9)/ sin(6'). 



which leads to the equivalent formula for the determinant 



det(a„(4)) 



n 



q — 2 cos n 



I 



i + 1 



;io.i9) 



l<l<i<n 

For instance, for n = 3, the determinant of the matrix Qslq) represented in Fig. 39 is equal 
to q^{q'^ — 1)'^((/^ — 2). The proof of the above formulas in Ref.[0 makes use of the intimate 
link between the meander problem and the Temperley-Lieb algebra ||6^. The equivalence 
between arch configurations and reduced elements of the Temperley-Lieb algebra TLn{q) 
goes as follows. The Temperley-Lieb algebra is expressed in its pictorial form, as acting 
on a "comb" of n strings, with the n generators 1, ei, 62, e^-i defined as 
1 



1 



1 

i+1 



(10.20) 

The most general element e of TLn{q) is obtained by composing the generators ( |10.20| ) like 
dominoes. The algebra is defined through the following relations between the generators 



(ii) [ei.ej] 
(Hi) Ci ei±i Ci 



qei z = 1, 2, n - 1 

if|z-j|>l 

e,- i = 1, 2, n — 1 



(10.21) 
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The relation (ii) expresses the locahty of the e's, namely that the e's commute whenever 
they involve distant strings. The relations (i) and (iii) read respectively 



O 



(iii) Cj Cj+i Cj 



1 

i+1 



(10.22) 

In (i), we have replaced a closed loop by a factor q. Therefore we can think of q as being 
a weight per circuit of string. In (iii), we have simply "pulled the string" number i + 2. 

An element e e TLn{q) is said to be reduced if all its strings have been pulled and all 
its loops removed, and if it is further normalized so as to read Yliei some minimal 

finite set of indices /. A reduced element is formed of exactly n strings. 





1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 



Fig.40: The transformation of a reduced element of TLg[q) into an arch con- 
figuration ^ of 9 arches. The reduced element reads = 636462656361666462. 



There is a bijection between the reduced elements of the Temperley-Lieb algebra 
TLn{q) and the arch configurations with n arches. Starting from a reduced element of 
TLn{q), we index the left ends of the n strings by l,2,...,n, and the right ends of the 
strings 2?i, 2n — 1, n-|- 1 from top to bottom (see Fig.40 for an illustration). Interpreting 
these ends as bridges, and placing them on a line, we obtain a planar pairing of bridges by 
means of non-intersecting strings, equivalent to a configuration of n arches. Conversely, we 
can deform the arches of any arch configuration so as to form a reduced element of TL^iq)- 
As a consequence, we have dim(TL^(q')) = c„, as vector space with a basis formed by all 
the reduced elements. 
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In the language of the Temperley-Lieb algebra, the meander matrix Qni^l) is precisely 
interpreted as the Gram matrix of the basis of reduced elements of TLn{q) (labeled by the 
corresponding arch configuration) with respect to the bilinear form (e^,ee) = Qn{(l)AB of 
Eq.( p.0.17|) . The determinant of Gn{(l) is obtained as a by-product of the Gram-Schmidt 
orthogonalization of this matrix, readily performed by use of the representation theory of 
the Temperley-Lieb algebra ||14|| . 

The formula ( |10.18|) may be extended to the more general case of meanders with open 
arches or meanders based on the Hecke algebra []7D[] . 

11. Two-flavor fully packed loops: exponents of the meander problem 
11.1. Generalized meanders as random lattice loop models 




Fig.41: An example of generalized meander with three loops of river (solid 
lines) and two circuits (dashed lines). The obtained graph is automatically 
bipartite. 

The most elaborate description of the meander problem consists in viewing it as the 
"gravitational" version (i.e. defined on random lattices) of a fully packed loop model. As 
in Section 10.3, we consider the general case of multi-river multi-circuit meanders (see 
Fig.41). We moreover attach a weight ni per loop of river and n2 per circuit. The case of 
multi-circuit meanders with a single river of Section 10.1 corresponds to taking ni — > 0, 
n2 = q, while that of multi-river multi-circuit meanders of Section 10.3 amounts to taking 
77-1 = 1 and 77-2 = q. 

We may therefore rephrase the problem as that of enumerating the configurations of 
planar graphs made of a set of river loops and circuits with the constraints that: 
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(a) The rivers are self- and mutually-avoiding; 

(b) The circuits are self- and mutually-avoiding; 

(c) A river and a circuit may cross at a "bridge" node; 

(d) The resulting graph is connected. 

To avoid problems of symmetry factors, it is as usual convenient to consider configura- 
tions with a marked and oriented edge (say corresponding to a portion of river). The above 
model is therefore a model of fully packed loops on random tetravalent graphs whose nodes 
correspond to the bridges, with two types of loops (the rivers and the circuits) weighted 
by ni and n2 respectively. We will denote this model by GFPL'^{ni,n2) (for gravitational 
fully packed loops). 

It is interesting to note that, as all nodes correspond to crossings of river loops and 
circuits, the corresponding tetravalent graphs are automatically bipartite, i.e. may have 
their nodes bicolored (see Fig.41). 



Fig.42: Examples of configurations of tangent meanders with two possible 
vertices, either the crossing of two loops, or a tangency point (contact with- 
out crossing). The underlying tetravalent graph may be arbitrary (a) or 
constrained to be bipartite (b). 

By anticipating the coming discussion, we may introduce two other versions of the 
problem, which consist in allowing for the presence, in addition to the crossings, of tangency 
points at which a circuit and a river come in contact but do not cross (see Fig.42). We 
shall refer to this case as tangent meanders. The bicolorability of the underlying graph 
is no longer ensured and, by analogy with the one-flavor loop case of Section 8, we may 




(a) 




79 



consider two different versions of the problem: the ordinary gravity version where we sum 
over aU tetravalent graphs or the Eulerian gravity version where we restrict the summation 
to node-bicolored tetravalent graphs only. 

11.2. The FPL"^ model on the square lattice 



X 



y 



Fig.43: Example of configuration of the FPL^ model. The black and white 
(here represented as dashed) loops may either cross or avoid each other. The 
two allowed vertices (up to rotations and symmetries) are represented on the 
right. 



We now would like to interpret the GFPL^{ni, n2) model or its generalizations with 
tangency points as gravitational versions of models defined on the regular lattice. In the 
case of tetravalent graphs, the regular lattice to be considered is the square lattice. A 
natural candidate is the -FPL^(ni, 77.2) model, introduced in Ref . [[fl| , which consists of a 
gas involving two types of loops, say black and white, with respective weights rii and n2 
and the constraints that: 

1. The loops are fully packed, i.e. each of the two systems of loops is self-avoiding and 
visits all the nodes of the lattice; 

2. Each link is occupied by one type of loop only, the loops being in contact only at the 
nodes. 

The FPL^(ni, n2) model allows for two types of vertices (up to rotations and symme- 
tries), represented in Fig.43 and corresponding respectively to a crossing of the two types 
of loops (crossing vertex) or to a contact with mutual avoidance (tangent vertex). Very 
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generally, we may introduce different weights, say x and y respectively for the two types 



of vertices. The FPL [ni, 72,2) model of Ref. concerns the particular case where x = y. 

Before we study this model, let us note that the suppression of the tangent vertex 
{y = 0) leads to a quite trivial model on the square lattice since the only allowed loop 
configuration is that where the black "loops" occupy the horizontal links of the lattice and 
the white loops the vertical links (or conversely). Still, by taking ni = 77-2 = 2, which 
amounts to orienting independently upwards or downwards each vertical (white) line and 
to the left or to the right each horizontal (black) line, the configurations that we obtain 
are in one-to-one correspondence with the foldings of the (dual) square lattice in d = 2, as 
studied in Section 1.1. The correspondence is simply that a vertical (resp. horizontal) line 
of the (dual) square lattice is folded if and only if the orientations of the two vertical (resp. 
horizontal) lines on each side are opposite. This remark, although elementary, prefigures 
the link between the GFPL^{2, 2) model and the folding of random quadrangulations. We 
will discuss this link more precisely in Section 13 below. 

Let us now consider the FPL?{ni, 77-2) model with x = y = 1. This model was studied 
in Ref. [|71| by Coulomb gas techniques similar to those presented in Section 5. Some 
properties could also be derived exactly by Bethe Ansatz fT^,^, in the same spirit as in 
Section 4. 

We start from the case ni = n2 = 2, which amounts to considering oriented loops 
with no extra weight. The model is then equivalent to a four-color model with colors, say 
A, B, C and D on the links of the square lattice and the constraint that the four colors are 
present around each node. The links of color A or B (resp. C or D) form the black (resp. 
white) loops and alternate along these loops. The two possible choices for the alternation 
of colors correspond to the two orientations of the loop. 

As in Section 3.3, we may transform the model into a model of heights X on the faces 
of the square lattice with the transition rules represented in Fig. 44. We need a priori four 
height differences A, B (black loops) and C, D (white loops), with the constraint that 
A ^ B + C ■\- D = 0to ensure a well defined height on each face. The obtained height 
variable X is therefore three-dimensional and corresponds to a three-component scalar 
field theory with central charge 3 JTlI . 



As in Section 5.2, the introduction of the weights ni and 77-2 is performed by intro- 
ducing vertex weights exp(±7rei/4) (resp. exp(±7re2/4)) per left or right turn of the black 
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Fig.44: Definition of tlie fieiglits X in tfie FPL^ model. Tlie definition uses 
tlie bipartite nature of the lattice. The consistency after one turn imposes 
the constraint A-\- B -\-C ^- D = 0. The height is therefore three-dimensional. 

(resp. white) loops, with ni = 2cos(7rei) and n2 = 2cos(7re2). For < ni,n2 < 2 
(0 < ei, 62 < 1/2), the theory is conformal with central charge ||71|| : 



CFPL(ni, n2) = 3 - 6 



6- 



1 - ei 1-62 



1 + Cdense 

(ni) + 

Cdense 



(11.1) 



with Cdense given by Eq.( |5.11|) . One can also compute the exponents xi^m associated with 
the long distance behavior r"^^' "' for the correlation of a set of / black lines and m white 



lines connecting two points at distance r, with the value ^1 

xi,m{ni, 712) 



1^^2_ 



2(1 -ei; 



;i-^^,o) 



,1-62 2 

H m — 



p2 
<^2 



2(1-62) 



{1-5. 



m,0) 



(11.2) 



-^Ia a (l-ei)(l-62) 

16 (1 - 61) + (1 - 62) 
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We may finally consider a "dense" version of the problem, denoted by DPL^{ni, 
(for densely packed loops) by allowing for the presence of sites not visited by black loops 
or white loops [|7^. In terms of height, this requires to impose the additional constraint 
A + B = = C + D, which leads to two-dimensional heights and reduces the central charge 
to a value 

CDPL(ni, n2) = 2 - 6- 6- — = Cdense(ni) + Cdense(n2) (11-3) 

1 - ei 1-62 

In this case, the exponents xi^rn reduce to 

, 1-62 2 (11.4) 
= xi{ni) + Xm{n2) 

with xi given by Eq.( |5.12| ). 

11.3. Coupling FPL'^{ni,n2) to gravity: Effective field theory of meanders 

The coupling of the FPL^(ni, 77-2) model to gravity corresponds stricto sensu to the 
case of multi-river multi-circuit tangent meanders with the two types of vertices, namely 
crossings and tangencies. For Eulerian gravity, realized by summing only over node- 
bicolored tetravalent graphs, the three-dimensional degrees of freedom of the regular lattice 
model are preserved and the critical behavior of the model corresponds to the coupling 
to gravity of a CFT of central charge cfpl (^1,^2)5 as given by Eq.( p.l.l|) . We may in 
particular use Eq.( [7.3|) to derive the configuration exponent governing the large n be- 
havior /i^ ~ p^"^ /n°' of the number /i„ of tangent meanders with one river line and one 
circuit meandering around it with 2n contact (crossing or tangency) points of alternating 
color along both circuit and river. This corresponds to taking ni,n2 —>■ 0, in which case 
cfpl(0,0) = —3, leading to 

a = ^^ (11.5) 



This is not however the original meander problem we started with. To get the 
GFPL'^{ni,n2) model, we have indeed to forbid tangency points. We shall now argue 
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X X+D-(A+D) 





X X+D-(A+D) 




(a) 



(b) 



Fig.45: Typical height configurations around (a) a crossing vertex and (b) 
a tangency vertex. When going diagonaUy say from the NW to the SE face 
(shaded), the height increases by an amount restricted to ±{A + C) or ±{A + 
D) in the case (a) and only to ih(A + B) in the case (b). Forbidding (b) 
will therefore amount to a reduction of the height variable range from 3 to 2 
dimensions. 

that tangency points are relevant and that their suppression reduces the value of the cen- 
tral charge. Assuming indeed that only crossing vertices are allowed, in which case the 
bicolorability of the nodes of the graph is automatic, we can still build a three-dimensional 
height as before. Note that the graph is also automatically face bicolorable (as is any 
tetravalent graph) and we may thus define a sub-lattice W of white faces and a sub-lattice 
B of black faces. It is then easy to see that the heights on two neighboring faces on the 
sub-lattice B, i.e. faces diametrically opposite around a node, may differ only by ±(A-|-C) 
or ±(A -|- D) (see Fig.45- (a)) but that the third direction ±(A -(- S) never occurs (as op- 
posed to what happens with the tangency vertex - see Fig.45-(b)). The heights of the 
sub-lattice B are de facto two-dimensional in the {A-\-C, A-\-D) plane and the same is true 
for those of the sub-lattice W . It is therefore harmless to set A + i? = C + D = 0, which 
takes us back to the DPL? model. This leads to the prediction that the multi-river 
multi-circuit meander problem lies in the universality class of the DPL'^{ni,n2) model, 
therefore a CFT with central charge 



Cmeander(ni,n2) = CdPl(?^1, ^^2) = Cdense(?^l) +Cdense(?^2) 



(11.6) 



The same conclusion may be reached by a different reasoning. Considering now tan- 
gent meanders on arbitrary tetravalent graphs (ordinary gravity version) , we may still use 
the height transition rules of Fig. 44 but we can no longer distinguish between A and —B 
(resp. C and — -D), which equivalently amounts to imposing A + B = Q = C + D too. 
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x+c A X x+c / X 




X A X-A X 



Fig.46: In the case of tangent meanders on arbitrary tetravalent graphs, the 
tangent vertices may be untied without affecting the heights. 

This also leads to a two-dimensional height and we again have to use the reduced central 
charge cdpl of Eq. (|11.6| ). To interpret this central charge as that for meanders, we still 
have to show that tangency points are irrelevant for ordinary gravity as opposed to the 
case of Eulerian gravity. We may understand this phenomenon heuristically by noticing 
that, in the effective height language, any tangent vertex may be safely "untied" as shown 
in Fig.46 without modifying the heights. 

To conclude, we have predicted in two different ways that meanders correspond to 
the coupling to gravity of a CFT with central charge given by Eq.( |ll.(^ ). For the original 
problem of meanders with one river line and one circuit, this gives c = —4 and leads to a 



configuration exponent |T6 



a ^ a(0) = 2 - 7.t.(-4) = ?^±^ (11.7) 



by use of Eq. (|7.3| 



11.4- More meander exponents 

We have now identified the CFT underlying meanders, as the dense two-fiavor loop 
model with 711 = ^2 = coupled to ordinary gravity. The complete knowledge of the 
conformal operator content of this CFT via the Coulomb gas picture gives access to a host 
of meandric numbers which we describe now. 

The important operators for our present purpose are those identified as creating ori- 
ented river vertices, namely the operators (pk (resp. (p-k), k = 1,2,... which correspond 
to the insertion of a /c-valent source (resp. sink) vertex at which k oriented river edges 
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originate (resp. terminate). In the Coulomb gas picture, these operators create "mag- 
netic" defect hues along which the height variable has discontinuities. The operator cj)^ 
has conformal dimension [ffll 



2 32 ^ ^ 

k = ±1,±2, with X as in Eq.( p.l.4| ). When coupled to gravity, these operators get 
dressed (into ^'fc) and acquire, according to Eq. ([7.4|) , the dimension: 



lv/8 + 3F-v/5 

29 -v^ ^ ^ ^ 



As a preliminary remark, we note that h±2 = A±2 = 0. The operators ^2-, ^-2 indeed 
correspond to the marking of an edge of the river in meanders, and moreover such operators 
must go by source/sink pairs for the orientations of the pieces of river connecting them 
to be compatible. Applying Eq.(|7.7|) to the two-point correlator {'^2'^ -2) A at fixed large 



area A = 2n, we find 

(*2^'-2)a ~ ^^f^ (11.10) 

while the meander counterpart (with a closed river) behaves as Ma/ (2A) ~ I (A^"'''^*''). 
We see that the net effect of the insertion of the operators \l/±2 is an overall factor propor- 
tional to oc n^, which confirms their interpretation as marking operators. 

We may now turn to the case of semi-meanders, for which the river is a semi-infinite 
line around the origin of which the road may freely wind. Considering the point at infin- 
ity on the river as just another point, the semi- meanders may equivalent ly be viewed as 
meanders whose river is made of a segment. Sending one of the ends of the segment to 
infinity just resolves the ambiguity of winding around either end. Using the above river 
insertion operators, we immediately identify the generating function for semi-meanders as 

(^'i^'_i) = J]5M,r7" (11.11) 

n>l 



Using again Eq.([7.7|) and the explicit values of A±i via (|11.9| ), we arrive at the large n 
asymptotics 



a ~" V 11 

^M^~^ a = l + 2Ai-7,t, = l + ^(v^+v^) (11.12) 
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star Eight | Cherry 

Fig.47: Three types of meandric configurations in which the river has the 
geometry of (a) a /c-valent star (b) an "eight" (c) a "cherry". The vertices 
corresponding to river sinks or sources are represented by filled circles (•). The 
edges of river in-between them are oriented accordingly. The road (dashed 
line) may freely wind around univalent vertices. 

Note that we expect the value of Qc to be the same for meanders and semi-meanders, as 
both objects occur as thermodynamic quantities in the same effective field theory. In the 
language of Section 10.1, this means that we are in regime (i). 

We may now generate many more meandric numbers by considering more general 
correlators. To name a few (all depicted in Fig.47), we may generate rivers with the 
geometry of a star with one /c-valent source vertex and k univalent sink vertices generated 
by (\I/fc(\E'_i)'^), rivers with the geometry of an "eight" with one tetravalent source vertex 
and two loops, each containing a bivalent sink vertex generated by (\1/4(\E'_2)^), or rivers 
with the geometry of a "cherry" with one trivalent source vertex, one univalent sink, and 
one loop, marked by a bivalent sink vertex generated by (\E'3\l/_i\l/_2), etc ... For each of 
these situations, we get the corresponding configuration exponent a = 3 — 7 -|- " 1) 

by applying Eg . (|777|) with the dimensions of Eq.( |11.9| ). We get respectively flf 



Q^k— star 



freight 



( V5 + y29)(V3/c2 + 8 + k{Vli - 2 V29) + 4^29 - 2^5) 
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^(v^-f v^)(v^-f v^) (11.13) 



ttcherry = ( + V^) ( + V^) 



11.5. Multi-circuit meander exponents 



More generally, we may consider the one-river multi-circuit case of Section 10.1 by 
taking tt-i ^ and n2 = q. Our predictions require moreover < < 2 for the field 
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theoretical description to hold. Using similar argument as in previous Section, we predict 
a critical exponent a{q) = 2 — 7str(c) with the central charge 

1 

c = Cmeander(0,g) = -1 + 6-^— , g = 2cos(7re), 0<e<- (11-14) 

J. c ^ 

We arrive at 



1 - e + 3e2 + a/(1 - e + 3e2)(13 - 13e + 3e^) . x . 

a(g) = 2 + A J g = 2cos7re 11.15 

d(1 — e) 

Similarly, we can compute the exponent o:{q) by inserting two defects corresponding 
to the two extremities of an open segment. This leads to the prediction a{q) = a{q) — 1 + 
2Ai(g) with Ai{q) related to hi{q) = xi,o(0,q)/2 of ([TL4D through (0). We finally get 

m 

. , , v/2(24e2 + e-l)(Vl-e + 3e"^+Vl3-13e + 3e^) wiiifi^ 
aw) = 1 + 24(1 -e) ^ " cos(7re) (11.16) 

The prediction for a{q) requires that q be less than a critical value given by = 2 cos(7rec) 
where Cc is the positive root of 24e^ + Cc — 1 = 0, namely [|T6 | 



,/g7 _ i\ 

gc = 2cos ( TT^-^^ — 1= 1.6738- •• (11-17) 

At this value, we have a —>■ 1. A heuristic argument shows that the entropically 
favored semi-meanders are then those with a "Russian doll" structure, i.e. those made of 
a first finite semi-meander using only the, say si first bridges closest to the source, of a 
second semi-meander using only the S2 next bridges (and possibly winding around the set 
source/first semi- meander), etc... with the s^'s all finite. We thus expect at qc a transition 
to a regime where the winding becomes extensive. The critical value qc is thus a good 
candidate for a transition point between the regimes (i) and (ii) of Section 10.1. 

12. Numerical checks 

The above predictions may be tested numerically by performing a direct enumeration 
of meanders, semi-meanders and other related configurations. Several algorithms have 
been used, the best ones producing results up to about 50 bridges. We shall present here 



two particular algorithms, one based on arch growth [|T5[, well adapted to the case of semi- 
meanders, and the other based on a transfer matrix formalism ||75| , [76| , p!7| . Beside exact 
enumerations, other statistical approaches such as Monte Carlo algorithms were also used 



in Ref. |77] to investigate large meander statistics. 
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transformation (I) 




Fig.48: The arch growth algorithm: the two transformations (I) and (II) are 
described in the text. X, Y et Z denote arbitrary systems of arches and Rn 
a rainbow made of n nested arches. Repeated application of (I) and (II) on 
the empty configuration generates all semi-meanders, represented here up to 
n = 4. 

12.1. Arch growth algorithm 

This algorithm generates configurations of semi-meanders. The meanders are then 
recovered as a particular case of semi-meanders with vanishing winding. We use the "open" 
representation of semi-meanders obtained by opening the river as in Fig. 38. A semi- 
meander with n bridges is therefore represented by a system of n arches above the (opened) 
river closed by a rainbow Rn made of n nested arches below the (opened) river. The 
algorithm generates all semi-meanders with n + 1 bridges from those with n bridges by 
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applying to the arch system one of the foUowing two transformations, as iUustrated in 
Fig.48: 

(I) we select one of the outer-most upper arches, cut it and re-close it below the river 
by encircling the original semi-meander, thus creating two arches above the river, two 
new crossings with the river and one new arch of rainbow below the river. In the 
original semi-meander picture, this transformation corresponds to the creation of one 
new bridge and does not modify the number of circuits. 
(II) we surround the existing semi-meander by a big circle, thus creating an outer arch 
above the river, two new crossings with the river and one new arch of rainbow below 
the river. In the original semi-meander picture, this transformation corresponds to 
the creation of one new bridge and increases by 1 the number of circuits. 
It is easy to see that the above algorithm generates all the semi-meanders with n + 1 
bridges once and only once from the set of semi- meanders with n bridges. Indeed, the 
transformation is easily inverted by cutting the external arch of the rainbow below the river 
and gluing it above. Each semi-meander is thus obtained from the empty configuration 
in a unique way by the application of a succession of transformations (I) and (II) (see 
Fig.48). An immediate consequence of this property is that the average number (out)„(q') 
of outer-most arches above the river for semi-meanders with n bridges weighted by a factor 
q per connected component is equal to 



Beyond mere enumeration, this algorithm has the advantage of allowing for a step 
by step tracking of the number of circuits, thus giving access to the analytic structure of 

(k) 

the SMn . More precisely, the number of circuits corresponds to the number of times 
the transformation (II) was used. For large n and k = n — I with / finite, we deduce for 
instance that 




R{q) - q 




The arch growth algorithm was used in Ref. |T^ to enumerate semi-meanders up to 27 



bridges. 



/ ! 



(12.2) 
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corresponding to the ('^) ~ n^//! choices of the transformation (I). More precisely, we have 

El 

^M^") = 1 
SM^-^^ = n - 1 n > 1 

SM'^-^^ = l(n2 + n - 8) n > 3 



(12.3) 



SM'^-^'^ = -in^ + Qn^ -?,ln-2A) n>5 
6 

^SM^"-^) = — (n'^ + 14n^-49n2 - 254n) n>7 

obtained by analyzing all possible shapes of semi-meanders with only 1, 2, 3, 4 applications 
of (I). 



12.2. Transfer matrix algorithm 




bridge #: 



operation: O U 



state: 



O D C* O C C 



q ^ c 
c 




Fig. 49: Transfer matrix algorithm. The operators O, C, U and D are de- 
scribed in the text. 



This particularly efficient algorithm was first used by Jensen in Ref. , where mean- 
ders were enumerated up to 2n = 48 bridges. The idea consists in generating the meanders 
from left to right one bridge after the other by means of a transfer matrix. The state of 



91 



the system between the m-th and the (m+ l)-th bridge is characterized by the connectiv- 
ity of those arches above and below the river which have been opened but not yet closed 
(see Fig. 49). These arches are indeed connected by pairs in the left part of the meander. 
These connections form themselves a transverse system of arches^ which, together with the 
relative position of the river, entirely characterize the state of the system. The application 
of the transfer matrix corresponds to the crossing of a bridge where one of the following 
four operations takes place (see Fig.49): 

O: Opening of a new transverse arch with extremities immediately on both sides of the 
river. 

C: Connection of the two extremities of transverse arches which are immediately on both 
sides of the river. If these two extremities belong to the same transverse arch, a new 
circuit of meander has been created on the left of that bridge. Otherwise, the number 
of circuits remains unchanged (operation denoted by C* in Fig.49). 

U: Upward migration (with respect to the river) of the system of transverse arches. 

D: Downward migration (with respect to the river) of the system of transverse arches. 

Starting from the empty state, the repeated application of the transfer matrix allows 
to generate meanders or semi-meanders according to the final state reached after 2n or n 
iterations. 



12.3. Numerical results 



The above algorithms lead to the following exact enumerations. The table below gives 
the numbers SM^^ of semi-meanders with n bridges and k connected components for n 
up to 27: 
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Semi-meanders with n bridges and k connected components 



174 

504 

1406 

4210 

12198 

37378 

111278 

346846 

1053874 

3328188 

10274466 

32786630 

102511418 

329903058 

1042277722 

3377919260 

10765024432 

35095839348 

112670458128 

369192702554 



140 
428 
1308 
4072 
12796 
404 32 
129432 
413900 
1342580 
4335288 
4201804 
46226896 
152594276 
500016036 
1660630740 
5472190206 
18264517264 
60475691308 
202684618564 
673892675030 
2266436498400 



125 

430 
1454 
4976 
16880 
57824 
197010 
675428 
2310268 
7927778 
27205180 
93448486 
321537086 
1105589516 
3812424912 
13121988240 
45330375774 
155172996170 
540314673578 
1853197292582 
5454265995454 



254 
956 
3584 
13256 
49052 
179552 
658560 
2394504 
8724464 
31575096 
114451388 
412811544 
1490190544 
5360943684 
19288139802 
69245171564 
243463024330 
890477645192 
3188033497580 
11409453277272 



104 

438 

1796 

7238 

28848 

113518 

444278 

1720384 

6643492 

25421620 

97136712 

368280210 

1395104236 

5250325378 

19746342212 

73863421394 

275113486146 

1027609557470 

3321478801772 

14161345139866 



3028 

12996 

54812 

228284 

333148 

3833076 



154( 



7428 



62244564 
247973928 
984221764 
3876113404 
15223550024 
59379645924 
231124139318 
394157177372 
3453279084296 
13266154255196 



1023 
4759 
21533 

415956 

1790C38 

7613504 

32094902 

134013518 

555915344 

2286993326 

9351284260 

33033821330 

153915109102 

518666292594 

2478777547126 

9370806527980 



51 
282 
1451 
7112 
33721 
155994 
708018 
3164400 
13365780 
60351384 
263582744 
1129603888 
4806762528 
20294554048 
35198877560 
355164255600 
1473618143854 
6074897243976 



62 
366 
1989 
10227 
50606 
243392 
1145250 
5295316 
24133016 
108659880 
484054494 
2136797594 
9351763592 
40532714144 
175243050328 
751285450888 
3199508532588 



74 
464 
2653 
14262 
73428 
366100 
1780690 
8492044 
39849468 
184478424 
344196480 
3823988336 
17168256704 
75446485224 
337955598100 
1433533303900 



87 

577 

3460 

19394 

103642 

534494 

2681612 

13162018 

63454876 

301374534 

1413224292 

5554314890 

30103132270 

137072634598 

519231327340 



70 6 
4428 



25320 

142940 

761120 

3932262 

19831956 

98069830 

477060460 

2238540120 

10347350938 

50377605328 

236416235332 



116 
852 



5576 

33758 

193274 

1061000 

5636950 

29174205 

147790613 

735416776 

3604466262 

17437690722 

83407238044 



132 

1016 

6924 

43448 

255380 

1451964 

7923525 

42037340 

217928136 

1108184612 

5543795384 

27346198448 



149 

1199 

8493 

55153 

336303 

1955009 

10947269 

59481467 

315291663 

1637050027 

8352021621 



16 
167 
1402 
10305 
59160 
434423 
2594686 
14895219 
82818768 
448523733 113659789 
2376167414 628492938 



17 

186 

1526 

12383 

85781 

554482 

3399516 

19990983 



1372 
14751 

105354 



34735627 



20 

249 

2434 

20458 

154844 

1084738 

7159268 



2752 
23850 
185576 
1333214 



3096 
27638 
220892 



Similarly, the table below gives the numbers MA ^ of meanders with 2n bridges and k 
connected components for n up to 20: 



Meanders with 2n bridges and k connected components 



42 
262 
1828 
13820 
110954 
933458 
8152860 
73424650 
678390115 
6405031050 
61505881612 
602188541928 
5959805669034 
59923200729046 
608188709574124 
6234277838531806 
64477712119534604 



640 
5236 
45164 
406012 
3772008 
35994134 
351173323 
3490581428 
35253449296 
360946635312 
3739935535756 
39159200588780 
413836299216608 
4409705753032543 
47337525317450815 
511563350415103008 



56 
580 
5894 
60312 
624240 
6540510 
69323910 
742513832 
8023001566 
87525544560 
961412790002 
10530964761765 
118257400015312 
1322564193693320 
14863191405246888 
167771227744292160 
1901345329556422790 



14 

240 

3344 

42840 

529104 

5413784 

76980830 

919032654 

10941339452 

130091632424 

1546154900644 

18379165280940 

218575190891815 

2601195918964184 

30980323101952368 

369239984556000356 

4405339231880790648 



42 
990 
17472 
271240 
3935238 
54787203 
742366152 
9871243396 
129477031190 
1531013354212 
21653861838390 
277215342450784 
3531450459003830 
44306855312849352 
556627894954179025 
7145814923879522986 



132 

4004 

85904 

1559934 

26200468 

412343728 

6230748192 

91335152243 

1310295995120 

18455675484196 

255301542207920 

3517938623373672 

47820311912336732 

6447533380737694 64 

8632733743310196255 



429 

16016 

405552 

8536390 

161172704 

2830421952 

47201994752 

757218376394 
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1860480 

44346456 
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5757601343920 
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1575533680321408 

25083534947894304 

390970503448090438 

5938061383039308348 
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5177542470 
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2175182185120 
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217139225832110534 

3587066097601934530 
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994840 
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27665550740 

637007339280 

13612373527632 

274491451519548 

5283147737822600 

97860442915614704 

1755310029771295216 
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143558599080 

3555368581553 

81730331353736 

1759714658172372 

36017597277512496 
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696717840 

24384586200 

727159530240 
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10850034322572432 
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2935292400 

112831907750 

3509538084152 

102530316254240 

2569370571766112 
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12713373760 

515709552000 
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532249561425750 
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942371200 
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84538447403968 
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10455495457248 
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951484123440 



1757263190 



705958959835306990 235265504762443572 64713641205591820 14658557362753320 2709804590253296 402058855155712 45500885566900 3978168315200 225760523500 6564120420 



These data may be analyzed in two different ways, either numerical, or analytic. The 
first approach, purely numerical, consists in building Mn{q) and SMn{q) for the values of n 
accessible from the data and to directly estimate R{q) (resp. R{q)) by taking for instance 
the ratio of two consecutive values of Mn{q) (resp. SMn{q)) and by using appropriate 
convergence algorithms to extract a limiting value at large n. 

The figure 50 shows the obtained estimates of R{q) and R{q) for values of q between 
and 6. One clearly distinguishes two regimes, a large q regime where R{q) > R{q) and a 
low q regime where R{q) = R{q) up to numerical errors (recall that ai q = 1, we already 
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Fig.50: Estimates of R{q) and R{q) for q between and 6. The numerical 
errors are smaller than the thickness of the curves. 

know that the two are indeed equal). This confirms the existence of a transition point 
(jc between a low q phase corresponding to the situation (i) of Section 10.1 of negligible 
winding, and a large q phase corresponding to the situation (ii) of extensive winding. The 
precise value of qc is difficult to estimate as the approach of the two curves is tangential. 
Still, the graph of Fig.50 is compatible with the analytic prediction of Eg. ( |1 1 . 17\ ) . 



a a V 




Fig.51: Estimates of the exponents a{q), a{q) and ^{q) for q between and 
8. The different curves correspond to different numbers of iterations of the 
convergence algorithm. 

We may similarly estimate the exponents a{q), a{q) and even u{q) by taking appro- 
priate combinations of (semi-) meander polynomials which converge to these exponents at 
large n. These estimates are represented in Fig.51 for q between and 8. The estimated 
value of u{q) has a large variation around q ^ q^ which suggest a discontinuity of u at 
the transition point. We therefore expect that u{q) varies continuously for q < qc (regime 
(i)) and has a jump discontinuity beyond which u{q) = 1 (regime (ii)). The exponent a(g) 
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follows the same scenario with a continuous variation for q < Qc and a jump discontinuity 
to a constant value a{q) = (we already know that a{q) 0). Finally, the exponent 
a{q) seems to decrease without discontinuity toward its limiting value a{q) 3/2. 

As far as the original problem of meanders and semi-meanders is concerned ^ 0), 



a refined analysis of the numerical data is presented in Refs.[^,|7§] and leads to the precise 
estimates 



R{q = 0) = R{q = 0) = 3.501837(3) 
a{q = 0) = 3.4208(6) a{q = 0) = 2.0537(2) 



(12.4) 
(12.5) 



These exponents agree remarkably with the predicted values a(0) = 3.420132 ... of 
Eq. dlLTl) and a = 2.053198 ... of Eq. dTTTl ). 




theory 




0.5 
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1 



(a) 



(b) 



Fig.52: Comparison between the theoretical value Cmeander(0, q) of the central 
charge and its estimated value from the numerical data of meanders (a) and 
semi-meanders (b). 



For generic q, we have represented in Fig.52 the comparison between the theoreti- 



cal central charge Cnieander(0, q') of Eq. (|11.14|) and its two numerical estimates: (a) that 
obtained by inverting the relation ( |11.15D from the measured value of a{q) and (b) that 
obtained by inverting the relation ( |11.16|) from the measured value of 6i{q). For g = 0, 
the estimates ( p.2.5| ) lead respectively in the cases (a) and (b) to c = —4.003(3) and 
c= -4.002(1). 

Note finally that the above transfer matrix method was also adapted to estimate con- 
figurational exponents for other meandric numbers such as those with "eight" or "cherry" 
river configurations (see Fig. 47). The agreement with the predictions (|11.13| ) is excellent 
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||17|| . Using the same method, it was also checked that the presence of the tangency ver- 
tex does not affect the exponents for the ordinary gravity version of meanders, hence it 
corresponds to an irrelevant perturbation, as we argued. 



A second way of exploiting the data consists in using them to obtain a large q expansion 
of various quantities. At large the expansion of R{q) up to order p in 1/q is governed 
by SM^ for / < p + 1. We have seen that for I finite and n large enough, SM^ is 
a polynomial of n of degree /. This property is a direct consequence of the arch growth 
algorithm. The existence of a thermodynamic limit for SMn{q)^^^ then guarantees that 
(y.{q) = as long as the 1/q expansion remains valid. Moreover, imposing that the above 
polynomials actually reproduce the numerical data leads to their complete determination 



up to / = 19. The net result is finally a large q expansion of R{q), given by [p!5 | 



2 2 2 4 8 12 10 4 12 46 

R{q) = 9 + 1 + - + — + ^- ^- -6 7~^~^ + ^ + ~ 

q q"^ q"^ q° q° q' q° q^ q^" q^^ 

98 154 124 10 102 20 64 ^( ^ ^ 

^12 ql3 gl4 ^15 ~ ql6 ^17 ~" ^18 ^^19 

A similar analysis for meanders allows to determine the analytic structure of the M, 
for large n and finite 1. Instead of polynomials, one obtains rational fractions in this case. 



;i2.6) 



(n-l) 



This results in fl5 



, ^/ 1 3 3 29 81 89 1 , , 

= 2V^(l + - + ^-^-^-^-^ + 0(^)) (12.7) 



7? 7? 




Fig.53: Comparison between the numerical estimates and the large q expan- 
sions for R{q) (a) and R{q) (b). For R{q), we have drawn the expansions 
truncated at order q~^ and q~^. For R{q), we have drawn the expansions 
truncated at order g"^, q~^'^ and q~^^. 
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The figure 53 shows the comparison between the values of R{q) and R{q) as estimated 
numericaUy and the expansions (|12.7D and (|12.6| ) above. 

The analysis of the corrections to the polynomial behavior for the numbers SM^ 
displays corrections of the form R{q)^ with R{q) given by Eq. (|12.70 . It is thus natural to 
think that the expansions above are valid as long as R{q) > R{q), i-e. in the whole regime 
(ii), and in particular that a{q) = 0, u{q) = 1 for all q > qc- For q < Qc, these expansions 
are no longer valid and therefore do not allow to describe the regime (i) . Fortunately, this 
regime (i) is precisely that for which we have obtained predictions from the KPZ formula. 



13. Meanders as folding of random quadrangulations 

In this paper, we started in part A by a study of the planar folding of the triangular 
lattice and extended it in part C to the case of random triangulations. We saw that 
this latter problem has a fully packed loop formulation on trivalent graphs. The above 
formulation of the meander problem as a fully packed loop gas on random tetravalent 
graphs suggests a link between this problem and the folding of quadrangulations. In 
this Section we will "close the loop" by showing that the meander problem is indeed the 
same problem as that of planar folding of random quadrangulations. This may be seen as 
the non-trivial gravitational version of the trivial correspondence between the FPL? {2^ 2) 
model with only the crossing vertex and the planar folding of the regular square lattice. 
This allows alternatively to view the folding of a self-avoiding one-dimensional chain as 
the phantom folding of a two-dimensional random quadrangulation, a somewhat surprising 
result. 

Let us start by considering the set of planar random quadrangulations, namely tessel- 
lations for which all the tiles are unit squares. By planar folding of such a quadrangulation, 
we again mean any map of the quadrangulation onto the plane whose restriction to each 
square tile is an isometry. The nodes of the quadrangulation then necessarily have their 
image on the nodes of a regular square lattice in the plane. 

o Foldability of random quadrangulations 

As in the case of triangulations, we first address the question of foldability in the 
plane. By use of the four-coloring of the nodes of the regular square lattice (see figure 
54), we immediately see that, if a folding exists, it naturally induces a four-coloring of the 
nodes of the quadrangulation. In this coloring, the four colors 1, 2, 3, 4 must appear in 
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Fig.54: Four-coloring of the nodes of the regular square lattice 




Fig.55: Example of a foldable quadrangulation (left), i.e. a quadr angulation 
whose nodes may be colored by four colors 1, 2, 3 and 4 in such a way 
that the nodes around each face are in cyclic order 1,2,3,4 or 4,3,2,1, also 
corresponding to an orientation of the links. The squares have been deformed 
in the planar representation. The equivalent dual tetravalent graph (right) is 
node-bicolored. The links dual to the 1-2 or 3-4 links form black lines which 
cross white lines formed by the links dual to the 2-3 and 4-1 ones. 

cyclic order 1,2,3,4 or 4,3,2,1 around each face. Conversely, a quadrangulation with such a 
four-coloring (which is unique up to a cyclic permutation of the colors) is clearly foldable 
onto a single square whose four vertices would have been labeled 1,2,3,4. We thus have 
the following equivalent characterizations for genus zero quadrangulations: 

1. The quadrangulation is foldable in the plane 

2. Its nodes may be colored by four colors 1, 2, 3 and 4 appearing in cyclic order 1,2,3,4 
or 4,3,2,1 around each face 

3. Its faces are bicolorable by two colors which must be distinct on two neighboring faces 

4. Its edges may be oriented in such a way that each face has a well-defined orientation 

5. The number of faces around each node is even 
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6. The number of edges adjacent to each node is even 

This last property 6 justifies that we call such quadr angulations Eulerian. 

The planar foldable quadrangulations are therefore the Eulerian quadrangulations, 
characterized by any of the above equivalent properties. 

The number of such quadrangulations with 4s edges (and rooted, i.e. with a distin- 
guished oriented link) is given by [[f8| 



.s-l 



AT, =8 



(2s + 2)!s! 60F si 



{3s) 



(13.1) 



to be compared with Eq. (|6.14| ). 



o Planar folding of random Quadrangulations 

The four-coloring of the nodes allows to classify all the links of the quadrangulation 
into the "horizontal" links 1-2 and 3-4 and the "vertical" links 2-3 and 4-1. These links 
are moreover naturally oriented in the direction 1— >2— s>3— s>4— >1 (see figure 55). The 
dual graph is a tetravalent node-bicolored graph, on which the dual of horizontal links 
form a gas of, say black, fully packed loops, and the dual of vertical links a complementary 
gas of white fully packed loops. The two systems of loops are further constrained to cross 
each other at the nodes of the tetravalent graph. 



T 



C 



.t 



1 
i 



=-t 



4 



D 

= +A 
= ±C 



Fig.56: Folding constraints for quadrangulations. In order for the folded 
image of a face to be a square of side unity, we must have ti = —ts = ±A 
and t2 — — ^4 = ±C where A and C are two orthogonal unit vectors. These 
constraints propagate to the neighboring squares so that each square has its 
two horizontal links equal to A and B = —A, and its two vertical links equal 
to C and D = -C. 
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It is then clear that, in a folded configuration of the quadrangulation, the faces are 
mapped on unit squares of the square lattice, hence the link variables (defined on the 
oriented links) may take only two values, say A and B = —A on the horizontal links, 
while those on the vertical links may take two values C and D — — C, where A and C are 
orthogonal unit vectors (see figure 56). Moreover, the two horizontal (resp. vertical) links 
around a given face have opposite values. 



11 10 




unfolded 



folded 



Fig. 5 7: Example of planar foldable quadrangulation (a) in a spherical rep- 
resentation. The corresponding black and white fully packed loops (b) are 
entirely fixed. Only their orientation remains arbitrary. The same loops in a 
planar representation (c) form a configuration of the GFPL^{2^ 2) model. The 
orientation of the white loops determines which horizontal links are folded, 
according to the rule displayed on the right. A similar rule determines which 
vertical links which are folded according to the orientation of the black loops 
(not represented here). For the particular configuration of loops drawn here, 
the links of the "perimeter" joining the nodes labeled 1 to 14 are necessarily 
folded irrespectively of the orientation of the black and white loops. The other 
links may be folded or not. For instance, the particular choice of orientation 
of the white loops represented here in (c) corresponds to having a horizontal 
fold between the nodes 12 and 8 in the representation (b), with thus four 
folded links. 



100 



In a dual representation, we can transfer the "color" A, B, C or D of the link variable 
to the dual link, in which case the colors alternate between A and B along the black loops, 
and between C and D along the white ones. Orienting each link of color A (resp. C) 
from the adjacent black node, and each link of color B (resp. D) from the adjacent white 
node, assigns a well-defined orientation to each loop. The orientations are independent on 
each loop, which precisely defines the GFPL^(2,2) model. The orientation of the loops 
corresponds to the folding degree of freedom, with the correspondence that a horizontal 
(resp. vertical) link is folded if and only if the two white (resp. black) loops immediately 
on each side of the dual link have opposite orientations (see figure 57). Conversely, the 
height variable associated to a configuration of the GFPL'^{2,2) model (which by defini- 
tion allows only for crossing vertices) with the rules of figure 44 and the choice B = —A 
and D = —C is nothing but the position of the nodes in the corresponding folded config- 
uration. To summarize, the GFPL^(2,2) model describes the planar folding of random 
quadr angulations. 

o Complete folding, partial folding and non-unfoldable quadran^ulations 

The situation is simpler than that of triangulations for the following reason. Given 
any foldable quadrangulation, or rather its dual tetravalent node-bicolored graph, there 
exists clearly a unique configuration of fully packed loops (up to a global interchange of 
black and white loops) compatible with the rule that loops have to cross each other at 
each vertex. For a given graph, there is no entropy associated with the positioning of the 
loops. The folding entropy is therefore entirely encoded in the choice of orientations for 
the loops. 

A first consequence is that the complete folding of quadrangulations, i.e. their folding 
onto a single square is simply described by the GFPL^(1, 1) model. Indeed, the folding 
onto a single square always exists and is unique for each foldable quadrangulation. On the 
tetravalent graph, it is obtained by the unique orientation of the loops defined as follows: 

1. We orient the outermost black loops counterclockwise; 

2. Inside each of these loops, we orient clockwise the outermost black loops; 

3. We continue iteratively for deeper loops by alternating the counterclockwise and clock- 
wise orientations; 

4. We repeat the procedure 1-3 for white loops. 

One clearly sees (see figure 58) that all the loops of a given color on neighboring nodes 
have opposite orientations, hence all links are folded. 
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Fig.58: Orientation of the loops leading to the complete folding. The black 
loops are oriented alternatively in the counterclockwise and clockwise direc- 
tion as we go deeper in the graph from the outermost loops. In this way, black 
loops on both sides of any white link always have opposite orientations, hence 
all vertical links are folded. The complete folding is obtained by orienting the 
white loops accordingly so that all horizontal links are folded as well. 

The foldable quadr angulations being in one-to-one correspondence with their complete 
foldings, we deduce that the GFPLP'il, 1) model also counts Eulerian quadrangulations, 
whose number is given by Eq. (|13.1D . 



Similarly, the GFPL^{1, 2) model describes the partial folding of quadrangulation, i.e. 
the complete folding in one direction only, with all vertical links folded, while the horizontal 
ones may be folded or not. The above interpretations are compatible with the values of the 
central charge Cmeander(2, 2) = 2 for the planar folding, Cmeanderll, 1) = for the complete 
folding, and Cmeander(l, 2) = 1 for the partial folding in one direction. 

The GFPL^{1, q) model was discussed in Section 10.3 with in particular the prediction 
( p.0.15|) - (|10.16|) for the asymptotics of the number of folded configurations. In particular 



for q = 1, one recovers i2(l,l) = 9/2as apparent on formula ( |13.1|) . For q = 2, the number 
of configurations of partial folding in one direction growths like (7r^/2)^* which, divided 
by the (9/2)^** foldable quadrangulations gives an average number of foldings varying as 
(tt^/Q)^'', i.e. (tt^/Q) per horizontal link. 

Finally, the original meander problem with ni = ?t.2 = may be interpreted as 
describing the foldable but "non-unfoldable" quadrangulations, i.e. the quadrangulations 
for which the complete folding is the only possible one. Indeed, requiring that no unfolding 
is possible amounts to requiring that all the relative orientations of the loops are fixed, 
which is true if and only if there is a unique loop of either species. 
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The mixed situation where ni = only amounts to selecting the quadrangulations 
which are non-unfoldable on vertical links, but possibly unfoldable on the horizontal links. 
For n2 = q = l, one imposes a complete folding of these horizontal links (although some 
of them could a priori be unfolded) while for 77-2 = = 2, one allows for unfolding the 
horizontal unfoldable links. 

To conclude, the passage via the meander picture has allowed to prove that the enu- 
meration of foldings on a single stamp of a closed, one-dimensional and self-avoiding strip 
of 2n stamps is equivalent to the enumeration of planar, two-dimensional phantom ran- 
dom quadrangulations with 2n faces, which are both foldable and "non-unfoldable", a 
quite amazing result. 
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CONCLUSION 



In this paper, we tried to give a unified presentation of a number of problems arising 
from the statistical physics of folding. We explained the intimate links between foldings 
of lattices or graphs, colorings of their links and/or nodes, as well as fully packed loop 
gases. Beside phantom folding, the introduction of self-avoidance in one dimension led us 
to the meander problem, itself related to phantom folding of random quadrangulations. 
We have introduced various formulations of these problems: constrained spin variables, 
height models. Coulomb gas description, all giving different angles of approach. 

In the case of regular lattices, the foldings in the presence of bending energy give rise 
to rich phase diagrams with conformational transitions intimately linked to the underlying 
geometry. In the case of random lattices, a similar sensitivity to the underlying graph was 
observed. It leads in particular to different universality classes according to the nature 
of the underlying lattice, distinguishing between ordinary and Eulerian gravity. More 
precisely, a shift of the central charge c — ^ c + 1 of the associated field theory is observed 
when going from dense loops to fully packed loops provided the latter are defined on 
bipartite lattices. 

This phenomenon seems to contradict the more familiar notion of universality in sta- 
tistical physics. Note however that a similar breakdown of universality occurs in other 
systems for which the symmetry of the problem is linked to the ambient geometry. A 
particular example is that of "hard objects" , say particles occupying the nodes of a lattice 
so that no two adjacent nodes be simultaneously occupied [179| , |80[| . On the regular lattice. 



these include the famous "hard square" problem on the square lattice , and the "hard 
hexagon" problem on the triangular lattice ||8^ (see also [Q). These two problems display 



a crystallization transition between an ordered phase at high density of particles and a dis- 
ordered phase of low density. The universality class of the transition point differs however 
for these two problems. This is easily understood by observing that the symmetry of the 
high density groundstates (corresponding to having particles on one sub-lattice made of 
next-nearest neighbors on the original lattice) differs. The square lattice is bipartite and 
has two such sub-lattices, leading to an Ising-like transition point, while the triangular 
lattice is tripartite with three such sub-lattices, leading to a transition point in the univer- 
sality class of the critical 3-state Potts model. When going to random lattices, a similar 



sensitivity was observed in the following sense: it was shown MM that hard particles on 



random tetravalent graphs display no crystallization transition when defined on arbitrary 
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graphs (ordinary gravity) while the transition is restored for bipartite graphs (Eulerian 
gravity). A simple heuristic explanation for this is that ordered groundstates exist only 
when the graph is bipartite. Remarkably, the ordinary or Eulerian nature of the graph is 
again what matters. 

Another surprise of our study is the appearance of irrational critical exponents for a 
priori simple combinatorial problems. Here again, other examples have been found. For 
instance, it was recently conjectured that the enumeration of closed planar curves with 
n intersections (a problem related to the enumeration of alternating knots) corresponds 
to the coupling to gravity of a conformal theory with central charge c = — 1, and that the 
asymptotic large n behavior is characterized by the same irrational exponent (|8.5|) as in 
Section 8. 

In these notes, we have limited ourselves to the simplest examples of folding. As we 
mentioned already, our description may be extended to the case of all compactly 2D foldable 
lattices of Fig. 5, or higher D-dimensional simplices, all with their own (multicolored) fully 
packed loop gas formulation and with their own symmetry. Similarly, meanders may be 
generalized in various ways: 



- meanders drawn on surfaces with higher genus [|T3[, accessible via topological expan- 
sion of the black and white matrix model of Section 10.3; 

- colored meanders, based on multicolor generalizations of the Temperley-Lieb algebra 

- SU{N) meanders, based on higher quotients of the Hecke algebra generalizing the 
Temperley-Lieb algebra (A^ = 2) [ffOj . 

Finally, a very natural question is that of a proper description of self-avoidance for 
lattices of dimension D > 2. At this time, only enumerations for rectangular domains of 



the square lattice with very small sizes are known |8^] and no systematic approach has 
been devised yet. 

The most advanced results presented here all rely on field-theoretical effective descrip- 
tions. It would be desirable to obtain the same results via a more traditional combinatorial 
approach. A first step in this direction was performed recently |^,57,?T| by solving various 
planar graph enumeration problems with combinatorics of trees. This approach is based 
on a proper cutting of the graphs into trees with local characterizations. Alternatively, 
it relates large graph statistics to properties of the so-called Integrated SuperBrownian 
Excursion (ISE), a probabilistic description of tree embeddings |8^. In view of these de- 
velopments, we may expect similar technique to apply to loop gas systems on random 
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graphs. Another probabihstic direction concerns stochastic Loewner evolution (SLE) in 
the context of which various exponents, predicted in Ref.[|8^ by use of the equation ( [7.3|) , 
were derived rigorously with the help of Brownian motion ||9^ . 

As a challenge for the reader, the simple combinatorial problem displayed in Fig. 31 
has a simple transfer matrix formulation and we may hope for an exact analytic solution 
to be within reach. 
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